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Abstract 

40(jm  by  40Mm  metal-barrier-superconductor  electron  tunneling  elements' 
responses  to  visible  laser  light  from  an  Argon  Ion  laser  are  examined. 
Through  an  in-depth  analysis  of  the  major  current  producing  processes,  it  is 
shown  that  once  the  impinging  E-field  is  coupled  across  the  barrier,  the 
dominant  response  results  from  the  modulation  of  the  Fermi  level  at  the 
applied  optical  frequency,  producing  a  rectified  current  due  to  the  nonlinear 
characteristics  of  the  inelastic  tunneling  process.  The  associated 
implications  of  this  with  respect  to  modulation  of  a  superconducting 
tunneling  current  at  an  optical  frequency  are  discussed,  leading  to 
applications  such  as  absolute  frequency  measurements  as  well  as  a  master- 
clock  oscillator  operating  at  an  optical  frequency. 
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Chapter  1 

1.1  Introduction 

Since  the  outset  and  formulation  of  quantum  mechanics,  many  of  the 
described  phenomena  clash  with  a  common  sense  (classical)  appreciation.  Of 
these,  tunneling  can  be  offered  as  the  most  striking.  Classically,  a  barrier  is 
a  barrier  not  to  be  trespassed:  a  ball,  with  an  insufficient  amount  of  energy 
to  scale  the  wall,  will  return  from  bouncing  off  the  wall.  The  trick,  of 
course,  is  to  use  really  tiny  little  balls,  and  lots  of  them.  Quantum  mechanics 
maintains  then,  that  if  the  wall  is  sufficiently  thin,  a  fraction  of  those  balls 
can  tunnel  through. 

Tunneling  followed  directly  from  Schrodinger's  equation,  and  was  used  to 
explain^  several  phenomena  that  were  coming  to  light  around  the  time  of 
Schrbdinger’s  original  paper,  such  as:  the  field  Ionization  of  atomic  hydrogen 
by  Oppenheimer  (1928),  field  emission  from  a  free-electron  metal  by  Fowler 
and  Nordheim  (1928),  alpha  decay  by  Gamow  (1928),  and  metal-insulator- 
metal  junctions  by  Sommerfeld  and  Bethe  (1933). 

Almost  thirty  years  later,  the  tunneling  phenomenon  fell  under  the 
scrutiny  of  the  engineers  when  in  1957  Esaki  developed  his  tunnel  diode, 
spurring  this  field  of  thought  to  a  new  height  of  interest.  Shortly  thereafter 
(1960),  Giaever  as  well  as  Josephson  studied  metal-insulator- 
superconductor  and  superconductor-insulator-superconductor  tunneling. 

Since  this  time,  many  experiments  have  been  accomplished  using  tunneling  as 
a  means  to  probe  the  superconductor's  energy  gap. 


'’For  an  exhaustive  history  and  references,  see  Solymar,  pp.  3.  19,  29-31,54-55. 
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More  recently,  however,  as  a  result  of  the  research  breakthroughs  at  the 
MIT  Optical  and  Infrared  Laser  Laboratory,  attentions  have  turned  toward 
high-speed  tunneling  junctions  with  applications  as  opticai  frequency  mixers 
[2][4-5][48-49].  These  metal-oxide-metal  point  contact  junctions  form  a 
class  of  high-speed  elements  with  the  novel  capabilities  of  responding  to  an 
applied  optical  field  at  the  applied  optical  frequency.  As  a  result  of  this 
work,  it  was  possibie  to  mix  microwave  and  far  infrared  radiation,  taking 
advantage  of  the  slight  nonlinearities  of  the  elastic  tunneling  (current 
flowed  through  the  contact  region  as  a  nonlinear  function  of  the  electric 
potential  applied  across  it).  The  first  optical  frequency  mixing  experiments 
were  performed  using  a  high  speed  adaptation  of  a  microwave  rectifying, 
metal-semiconductor  (forming  a  Schottky  barrier),  point  contact  diode[4], 
successfully  generating  harmonics  of  microwave  radiation  mixing  with  337 
pm  (-10^^  Hz)radiation.  This  could  occur  because  the  capacitor-like  junction 
could  oscillate  at  these  frequencies. 

However,  the  intrinsic  time-constant  limitations  of  the  semiconductor 
element  necessitated  a  different  approach  for  higher  speed  work.  This 
resulted  in  poiished  metal  posts  replacing  the  semiconductor  [5],  forming  a 
metal-barrier-metal  junction,  which  were  then  responsible  for  frequency 
mixing  out  to  1.5  pm[48].  But  the  associated  inherent  vibration  sensitivity 
of  this  geometry  precluded  its  further  study  in  the  search  for  faster 
responses. 

In  order  to  better  understand  and  to  study  the  physical  processes  of 
electron  tunneling  of  the  metal-barrier-metal  junctions  lead  to  the  geometry 
of  the  thin-film  metal-crosses  that  were  used  by  Elchinger  [1]  and  were  used 
in  this  work,  and  whose  dimensions  were  on  the  order  of  tens  of  microns. 
These  crosses  are,  in  effect,  a  tiny  capacitor  formed  by  evaporating  a  thin. 
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metal  strip  onto  a  quartz  substrate,  oxidizing  it,  and  then  sandwiching  this 
oxide  layer  with  the  second  metal  strip.  Judiciously  choosing  one  of  these 
metals  such  that  it  goes  superconducting  at  a  convenient  cryogenic 
temperature,  it  is  then  possible  to  take  advantage  of  the  inherent  and  much 
more  significant,  tunneling  nonlinearities  due  to  the  band-gap  effect  of  the 
metal-oxide-superconductor  electron  tunneling  element. 

Until  the  present  work,  several  investigators  have  studied  the  response 
of  these  same  type  junctions  to  electric  fields  whose  frequencies  were  at 
about  the  (RC)-"'  time  constant  of  the  junction  as  well  as  at  higher 
frequencies  [6-10].  They  concluded  that  the  lower  frequency  (in  the  Rf) 
response  mechanism  arose  from  the  nonlinearities  of  the  junction's  l-V 
characteristics  and  was  due  to  rectification  as  a  result  of  the  modulation  of 
the  Fermi  level  across  the  barrier,  but  that  the  higher  frequency  (in  the  far 
and  near  IR)  response  was  due  to  a  photo  induced  tunneling  [1]. 

As  a  result  of  this  work,  It  can  be  shown  that  the  junction's  response  to 
the  higher  frequency  (up  to  the  IR  and  visible),  contrary  to  the  earlier  works, 
is  indeed  due  to  modulation  of  the  Fermi  level  across  the  barrier  at  the 
applied  frequency,  causing  rectification,  and  that  this  process  can  be  made  to 
dominate  over  others.  This  is  important  in  that  it  is  possible  to  take 
advantage  of  the  higher  order  of  nonlinearity  (those  far  exceeding  that  of  a 
normal  metal-barrier-normal  metal  junction)  in  order  to  accomplish 
frequency  mixing  well  into  the  visible  and  even  into  the  UV. 

1.2  Thesis  organization 

The  purpose  of  this  thesis  is  to  examine  the  response  of  a  metal-barrier- 
superconductor  when  it  is  exposed  to  light.  Several  processes  are  involved, 
some  of  which  are  competing,  while  some  have  nearly  the  same  response 
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shape.  The  matter  at  hand,  then,  is  to  show  under  what  conditions  each 
process  dominates. 

There  are  several  steps  in  making  this  identification,  and  as  such,  lend 
themselves  nicely  to  the  thesis  organization; 

•  Chapter  2  outlines  the  basic  theory  of  tunneling,  to  include  the 
1-dimensional  basics,  metal-insulator-metal  tunneling,  and  metal-barrier- 
superconductor  tunneling.  In  addition,  the  processes  of  photoemission, 
rectification  resulting  from  Fermi  level  modulation  across  the  barrier,  and 
how  a  differential  temperature  change  (due  to  the  junction  being  heated)  can 
manifest  itself  in  an  observed  signal  are  all  discussed. 

•  Chapter  3  describes  the  experimental  setup  and  the  procedures  used  to 
collect  the  data. 

•  Chapter  4  details  the  data  analysis  in  determining  the  associated 
barrier  parameters,  compares  theoretical  predictions  with  experimental 
results  for  the  l-V  curves  as  well  as  the  associated  derivative  curves,  and 
accommodates  earlier  interpretations. 

•  Chapter  5  deals  primarily  with  the  implications  of  this  work, 
improvements,  and  future  experiments  that  can  and  will  dovetail  from  these 
findings  to  include  more  efficient  ways  to  couple  the  light  across  the 
junction,  how  to  make  higher  speed  junctions,  and  the  future  of  frequency 
mixing  and  absolute  frequency  measurements. 

•  Finally,  several  appendices  have  been  included  for  ready  reference  and 
clarity.  The  first  appendix  details  from  a  microscopic  view  the  constancy  of 
the  superconductive  calculations  with  the  corrective  factor  of  a  change  in 
density  of  states;  the  second  appendix  outlines  the  advantage  of  four-pole 
impedance  measurements,  and  the  third  appendix  describes  the  Cray  X-MP 
FORTRAN  subroutine  used  to  numerically  evaluate  the  l-V  integrals. 
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Chapter  2 

Tunneling  Theory 
2.1  One-Dimensional2 

Much  insight  and  relevance  can  be  gleaned  about  the  processes  of  interest 
in  this  experiment  by  starting  from  basic  principles,  and  examining  even  the 
one-dimensional  case. 

An  electron  moving  with  a  kinetic  energy  of  E  in  zero  potential  is  incident 
upon  a  potential  barrier  U2,  such  that  U2>E,  as  shown  in  Figure  1. 


I 

II 

III 

E 

U2 

- > 

U3 

C 

II 

0 

1 

1 

1 

0  w  X 

Figure  1  —  One-dimensional  potential  barrier. 


Now,  from  the  time-independent  Schrodinger  equation,  we  car.  work  out  the 
probability  of  the  electron  appearing  on  the  other  side  of  the  barrier 


__B_V^  +  uU  =  Eip 

2m 


(2.1.1) 


^theory  for  this  chapter  has  been  taken  from  the  following  references;  [11],  [12],  [3,  pp  20-21). 
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.  or  a  constant  potential  in  region  I,  we  can  find  the  wave  equation  of  the 
electron  to  be 


where 


tp,  =  Aexp(/kjx)  +  Bexp(-/k]X) 


(2.1.2) 


(2.1.3) 


Now  the  solutions  in  regions  II  and  III  are 


where 


H>2  =  Cexp(-k2x)  +  Dexp(k2x) 
=  Fexp(/k^x) 


k2  =  ^;2m(U  2  -  E  ^  and  k^  =  1  2m(U  3  -  E  ) 


(2.1.4) 

(2.1.5) 


(2.1.6) 


and  A,  B,  C,  D,  F  are  constants  which  we  can  determine  by  imposing  the 
condition  that  the  wave  function  is  smooth  and  continuous.  Hence,  at  the 
boundaries  x  =  x,,  and  x  =  Xj,  ^  and  dH’/dx  must  be  continuous.  Now  the 

quantity  of  real  interest  is  |F/A|,  the  amplitude  relationship  between  the 
output  and  the  input.  Assuming  that  k2W  »  1,  then  exp(-k2W)  is  negligible 
compared  to  exp(k2w),  and  it  is  possible  to  show  that 


4k, k2 


k,  +  k2pjk3  +  k2 


exp(-k2w) 


(2.1.7) 


Now,  we  can  obtain  the  ratio  of  current  densities,  such  that 


),  “lA  (k^.k^jjk^k^) 


(2.1.8) 


What  is  of  extreme  importance  here  is  that  the  above  expression  is 
dominated  by  the  exponential  factor.  Consequently,  the  wider  the  barrier  is, 
the  smaller  the  current  that  can  flow  across  it.  From  many  previous  studies, 
practical  barriers  in  solids  tend  to  be  on  the  order  of  1  eV.  Using  this  value 
for  (Uj  -  E)  in  Equation  (2.1.6),  taking  the  free  electron  mass,  and  a  current 

transmission  of  10‘®,  then  the  width  of  the  barrier  should  be  on  the  order  of 

w  =  InlO^- lO'^m  (2.1.9) 

2k2 


in  order  to  observe  tunneling. 

2.2  Metal-Barrier-Metal  (MBM)  Tunneling^ 

Before  discussing  metal-barrier-superconductor  tunneling,  it  would  be 
best  to  lay  the  groundwork  with  first  a  short  investigation  of  metal-barrier- 
metal  tunneling.  From  a  simple  model  (Figure  2a),  it  is  possible  to  see  that 
in  thermal  equilibrium  the  Fermi  levels  of  the  two  metals  are  equal.  Now, 
defining  a  net  positive  electron  flow  as  being  from  left  to  right,  when  a 
negative  voltage  is  applied  to  the  metal  on  the  left,  all  the  electrons  of  that 
metal  will  have,  then,  an  energy  eV  higher.  As  a  result,  the  relative  Fermi 
energies  can  be  depicted  as  in  Figure  2b. 


^theory  for  this  chapter  has  been  taken  from  the  following  references:  [131.(14],  (15, pp  31-46], 
(3,  pp.  22-23], 
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Figure  2  -  Fermi  energy  diagram  for  a  metal-barrier-metal 


junction. (a)  Thermal  equilibrium;  (b)  at  a  potential  difference  of  V. 

There  are  several  approaches  to  calculate  the  current  flowing  across  the 
junction,  but  here  it  will  be  sufficient  to  say  that  the  number  of  electrons 
moving  from  left  to  right  in  an  energy  interval  dE  must  be  proportional  to  the 
number  of  occupied  states  on  the  left.  In  other  words, 

N,(E  -eV)f(E  -eV)dE  (2.2.1) 

where  N  i  is  the  density  of  states  of  the  metal  on  the  left,  and  f  is  the  Fermi 
function  (taking  the  reference  point  as  being  the  Fermi  level  of  the  metal  on 
the  right,  since  in  the  experiment,  it  is  the  grounded  metal).  Now,  for  there 
to  be  a  current  moving  from  left  to  right  there  must  be  unoccupied  states  in 
the  metal  on  the  right.  In  other  words,  the  current  must  also  be  proportional 
to 

N^E)[l-f(E  -eV)]dE  (2.2.2) 
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where  is  the  density  of  states  of  the  metal  on  the  right.  Incorporating  all 
the  above  with  the  probability  of  transition  across  the  barrier,  Pi-t(E)  ,  the 
current  flowing  from  left  to  right  is 

I^„-Pl^iE)N,(E  -eV)f(E  - eV)N,(E )[l-f(E  -eVjdE  (2.2.3) 

while  the  current  flowing  from  right  to  left  is 

I,,|-P.-(E)N,(E  -eV)f(E  -eV)N,(E)[l-f(E  -eV)jdE  (2.2.4) 

Making  the  assumption  that  P|-»r{E)=  Pr-H(E) ,  (that  is,  an  electron  has  just  as 
much  chance  to  tunnel  from  left  to  right  as  from  right  to  left),  and 
integrating  over  all  energies,  the  net  current  becomes 

I'j  P,^JE)N,  (E-eV)N,(E)(f(E  -  eV)-f(E)l  dE  (2.2.5) 

Now,  interestingly,  for  small  applied  voltages,  it  is  not  uncommon  to  take 
Pi-.r{E}  as  being  independent  of  energy,  thus  removing  it  from  under  the 
integral.  Similarly,  it  can  be  argued  that  the  density  of  states  is  a  slowly 
varying  function,  and  simply  take  its  value  at  the  Fermi  level,  or 

N,(E-eV)*Ni(E)=N,(0)  and  N^E)«N^0)  (2.2.6) 

Incorporating  the  geometry  considerations  of  the  barrier/junction  as  well  as 
the  PHr(E)  into  A,  the  expression  for  the  current  reduces  to 

I  .  A  N,(0)N,(0)[ [f(E-  eV)-  f(E  )]dE  (2.2.7) 
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This  then  allows  the  Fermi  functions  to  be  expanded  for  small  voltages  to 
yield 


[f(E  -eV)-f(E)j=  (2.2.8) 

dE 

which,  if  the  temperature  is  not  too  high,  can  be  approximated  by  a  delta 

dE 

function,  giving 

1  =  AN,(0)Nj0)eV  (2.2.9) 

In  short  then,  a  metal-barrier-metal  junction  obeys  Ohm’s  Law,  with  a  linear 
relationship  between  current  and  voltage. 

The  above  heuristic  approach  was  convenient,  but  to  validate  its 
simplicity  a  basic  set  of  states,  a  perturbing  Hamiltonian,  a  final  set  of 
states,  and  a  matrix  element  connecting  the  two  states  should  be  analyzed. 
This  approach  was  actually  tackled  by  Bardeen  [16]  first,  then  later  by  Kane 
[15,  pp.  1-12]  in  greater  detail. 

^2  in  Equation  (2.1.4)  represents  an  exact  solution  of  Schrodinger's 
equation  in  region  ll,  where  there  are  two  exponential  functions  decreasing  in 
the  positive  and  negative  directions.  Rewriting  these  two  functions 

ip,  =  C),exp (-kjx)  and  il)^=C^exp(-k2x)  (2.2.10) 

Now,  assume  that  is  a  solution  not  only  for  x^  <  x  <  X2,  but  for  the  whole 
region  x  >  x^.  This  incurs  only  a  slight  error  since  is  already  extremely 
small  at  X  =  X2.  Similarly,  consider  as  a  solution  for  x  <  X2.  Assuming, 
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then,  that  an  electron  is  in  the  initial  state  of  it  is  possible  to  work  out 
its  transition  rate  into  complete  time-dependent  wave  function 


^  =  q(t)  %e\p  |-*  +  ’tt)  I"' (2.2.11) 

Substituting  the  above  equation  Into  the  time-dependent  Schrodinger 
equation 

(2.2.12) 

dt 

and  assuming  that  only  a  very  small  proportion  of  electrons  make  the 
transition  (thus,  implying  q  =  1,  and  r  =  0),  and  that  the  electrons  are  either 
in  'P,  or  ,  it  follows  that  dq/dt  «  0.  Using  this,  multiplying  both  sides  of 

the  equation  by  'P*  and  integrating  yields 

ifi^  =  J  -  E  )q>,dxexp|/5tl^|  (2.2.13) 

Writing  the  Hamiltonian  H  in  the  form  of  Ho  +  Hi  and  noting  that  Hq'Pi  =Ei'Pi 
yields 


(2.2.14) 


SO  that  now  the  effective  tunneling  matrix  element  may  be  defined  by 


I 


-  E,)ip]dx 


(2.2.15) 
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Now,  in  review,  the  essentiai  assumptions  in  the  above  derivation  were  that 
and  were  good  approximate  solutions  of  the  exact  Hamiltonian  H,  and 

exact  solutions  of  the  hypothetical  'unperturbed'  Hamiltonian  Hq,  and  thus 
H'P,  =  E'l’i  for  X  <  x2,  and  H4’^  =  for  x  >  x1 .  Hence  the  integral 


Xb>Xi 


(2.2.16) 


is  zero,  and  may  be  added  to  the  integral  in  Equation  (2  2.15)  yielding  the 
symmetric  form 


-El  jipi-iiJilH-E^jip^jdx 


Integrating  by  parts  yields 


XiSXbSX2  (2.2.17) 


(2.2.18) 


which,  upon  substituting  Equations  (2.2.10)  into  the  above  gives 

Trt  =  ^c;c,  (2.2.19) 

With  this  matrix  element,  the  current  may  be  determined  by  using  Fermi's 
Golden  Rule,  yielding  the  same  expression  that  was  obtained  for  the 
transmitted  current  before.  This  is  critical,  for  it  shows  that  the  Bardeen 
approach  provides  the  same  result  as  the  simple  one-electron  calculation. 
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2.3  Metal-Barrier-Superconductor^ 

Using  the  usual  semi-conductor  representation,  Figure  3  shows  a  junction 
formed  by  sandwiching  a  thin  insulating  barrier  between  a  normal  metal  and 
a  superconductor.  At  absolute  zero  temperature,  ait  of  the  states  in  the 
superconductor  are  filled  up  to  Ep  -  A,  with  no  filled  states  above  the  gap. 


(a)  (b) 

Figure  3  --  Energy  diagram  of  a  metal-barrier-superconductor  junction 

at  T=  O’K,  (a)  V  =  0;  (b)  V  >  A/e 


When  in  thermal  equilibrium,  the  Fermi  energies  must  match  (Figure  3a). 
When  a  voltage  V  <  A/e  is  applied,  the  electrons  on  the  left  have  no  empty 
states  on  the  right  to  access,  and  thus  no  current  flows.  At 
V  =  A/e,  there  is  a  sudden  increase  in  current,  not  only  because  now 
electrons  may  tunnel  from  left  to  right,  but  they  also  face  a  large  density  of 
available  states.  As  the  voltage  is  Increased  further,  even  more  empty 
states  become  available  for  tunneling  and  the  current  increases  as  indicated 
in  Figure  4.5 


^theory  for  this  chapter  has  been  taken  from  the  following  references:  [3,  pp.32-391,  [17],  [18]. 
5 For  an  extremely  palatable  account  of  the  first  works  in  this  area,  see  [17], 
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Figure  4  -  The  l-V  response  from  a  metal-barrier-superconductor 

at  T=0*K 

Now,  at  finite  temperatures  (Figure  5),  there  are  electrons  above  Ep  +  A, 

as  well  as  holes  below  this  gap.  In  addition,  some  of  the  electrons  of  the 
normal  metal  are  in  excess  of  Ep  +  A,  even  at  thermal  equilibrium. 

Consequently,  even  a  tiny  voltage  Is  sufficient  to  initiate  a  current,  but  a 
sizable  rise  in  current  will  again  occur  at  V  =  A/e  (Figure  6) 


Figure  5  -  Energy  diagram  for  a  metal-barrier-superconductor  at  T>0“K, 

in  thermal  equilibrium 
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With  this  background,  and  the  formaiism  of  section  2.2,  it  is  possible  to 
develop  an  expression  for  the  current  that  will  flow  through  a  metal-barrier- 
superconductor  junction.  The  essential  assumption  is  that  when  the  metal 
goes  superconducting,  only  the  density  of  states  changes,  but  nothing  else 
does  (the  matrix  element  remains  the  same)®.  Consequently,  Equation 
(2.2.7)  may  be  used  to  describe  this  current,  provided  the  density  of  states  of 
the  superconductor  is  taken  into  account.  That  is,  replacing  Nr(E)  by  the 
density  of  states  in  the  superconductor,  Ng, 

Ns(E)  =  N,^E)n^E|  (2.3.1) 

where  is  the  density  of  states  of  that  metal  in  the  normal  state,  and  is 
the  BCS  superconductor,  energy  gap-dependent,  density  of  states  correction 
term.  Now  the  current  may  be  written  as 

I„s-AN„,(0)N^0)j  «,|E)'f(E-eV)-f(E)'idE  (2.3.2) 

where  Ni^andN^jsi  refer  to  the  normal  metal  density  of  states  of  the  metals 
on  the  left  and  right,  respectively,  in  Figure  5.  At  T  =  0,  the  Fermi  function 
is  unity  up  to  f(0),  and  takes  the  value  zero  above  that.  This  then  modifies 
the  above  expression  for  the  current  to 


®See  appendix  1. 
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/•eV 


Ins  =  ANJO)N^O) 


J  A 


EdE 


E^-A^ 


1/2 


(2.3.3) 


=  AN^O)N^O)[(eVf 


(2.3.4) 


for  eV  >  A,  while  0  for  eV  <  A. 

Now,  when  T>0,  it  is  necessary  to  integrate  over  all  possible  energies 
so  that 


Ins=ANJ0)N^0) 


|E|[f(E-eV)-f(E)jdE 


E^-A^ 


1/2 


(2.3.5) 


Although  this  has  been  solved  in  closed  form  for  eV  <  A  (with  some 
simplifying  assumptions)  by  Giaever  and  Megerle  [19],  the  predominant 
region  of  this  present  work's  interest  lies  just  prior  to  where  the  junction's 
response  goes  Ohmic.  As  a  result,  numerical  integrations^  were  required  to 
calculate  the  predicted  curves  for  the  metal-barrier-superconductor 
junctions  that  were  used. 


2.4  Current  Producing  Responses^ 

2.4.1  Photoemission 

Now  that  the  ground  work  for  the  experimental  conditions  of  the 
tunneling  element  has  been  set,  this  leads  to  the  discussions  involving  the 
mechanisms  that  yield  current  producing  responses.  From  the  vast  amount  of 

^see  appendix  #3. 

^theory  for  this  section  has  been  taken  from  the  foiiowing  references;  [1],  [14],  [29]. 
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work  that  has  already  been  done  in  this  area  [20-28]^  the  dominant  current 
producing  junction  responses  are  due  to  photoemission,  rectification,  and  a 
thermally  induced  current. 

Photoemission  and  its  bevy  of  intricacies  have  been  studied  extensively. 
Applying  the  formulations  to  this  tunneling  junction  environment,  Stratton 
[14],  Burshtein  and  Levinson  [29],  Elchinger  [1],  and  Jaklevic  and  Lambe  [30] 
have  already  compiled  the  detailed  derivations.  Following  predominantly 
Stratton's  approach,  the  equation  describing  the  photoemissive  current  for  a 
trapezoidal  barrier  (Figure  7)  can  be  obtained. 


Figure  7  —  Trapezoidal  Barrier 
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Amending  the  one-dimensional  approach  to  contain  electrons  moving  in 
all  three  axial  directions,  the  current  density  can  once  again  be  written 
phenomenologically  as  the  product  of  the  flux  of  normally  incident  electrons 
on  the  barrier  at  energy  Ex  potentially  contributing  to  the  current,  and  the 
probability  of  penetration  through  the  barrier,  integrated  over  all  applicable 


J  = 


^-eVu 


F(E,,VjP(E,,V,)dE, 


(2.4.1) 


Simmons  computes  the  electron  current  density  per  unit  energy  (the  supply 
function)  by  averaging  the  electron  energy  density  parallel  to  the  barrier 
(i.e..  In  the  y  and  z  directions),  leaving  only  the  x-directed  flux.  Consequently, 

r 


F(Ej  =  4l5mfi  f(E„E,)dE, 

h  Jo 


where 


f(E«.E,)  = 


J 


1 


(2.4.2) 


(2.4.3) 


Now,  integrating  over  all  radial  energy,  Er, 


(2.4.4) 

h 

where  e  is  the  electronic  charge,  h  is  Planck's  constant,  k  is  Boltzmann's 
constant,  T  the  temperature  of  the  electrode,  and  &  is  merely  1/kT.  Just  as 
before,  the  net  current  will  simply  be  determined  by  the  difference  in  supply 
functions  induced  by  an  external  bias  voltage,  Vb.  such  that 


j  j  _4imie(kT) 

•'Tot  •’l— r  •'i-.l  i 


dE. 


(2.4.5) 
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Now,  unlike  the  rectification  response  where  the  Fermi  level  at  the 
barrier  oscillates  at  the  applied  electric  field  frequency,  the  photoemissive 
current  is  a  DC  effect,  where  only  a  fraction  of  the  electrons  in  the  Fermi 
sea  contribute.  Assuming  that  the  absorbed  light  will  be  equally  partitioned 
among  the  three  orthogonal  velocity  components,  and  that  initially  the  light 
illuminates  only  the  normal  metal  side  (the  left  side  in  Figure  7),  ther.  at  Vb 
=  0  only  electrons  with  an  energy  'hv  will  contribute  to  the  current.  This 
results  in  having  to  shift  the  supply  function  from  Ex  to  (Ex  -  hv)  in  order  to 
account  for  these  electrons.  Of  all  the  electrons  that  are  photoexcited,  only 
a  fraction,  f,  of  them  will  have  their  velocities  directed  in  the  x-direction. 
As  a  result  then,  the  photocurrent  can  be  written  as 


4jime  (kT) 


'P.E. 


1— r 


/ 


f  P(Ejta( 

Jo 


1  + 


(2.4.6) 


During  this  experiment  P-  1 /!(0.086  meV/‘*K)(2®K)  »  1  eV"’  while  the 

photon  energies  were  2.7  eV  or  less,  allowing  the  approximation 

/ 


|0. 

_4!!m£(E  -hv), 
h’ 


for  hv 
for  <  hv 


(2.4.7) 


As  will  be  shown,  hv  is  above  the  barrier  height  in  all  frequencies,  and 
assuming  that  attenuation  by  scattering  or  through  other  energy  loss 
schemes  is  negligible,  it  then  becomes  possible  to  assume  P(Ex)  «=  1. 

Allowing  for  the  superconducting  side  to  contribute  a  fraction,  g,  of  the 
electrons  to  the  net  current,  the  total  photocurrent  becomes 


25 


(Ex-HdE. 


(2.4.8) 


J 


P.E. 


“l(/‘-g)(hv-cp2f 


(2.4.9) 


where,  ^ 

It  is  from  this  last  expression,  that  the  barrier  height,  ‘P2,  is  determined 
from  a  Fowler  plot  [31]  by  plotting  versus  hv,  (see  Figure  17),  where  ^2 

is  the  x-axis  intercept 

2.4.2  Rectification 

Now,  another  conduction  process  characterizing  the  fundamental  response 
of  the  metal-barrier-metal  junction  to  electromagnetic  fields  is 
rectification,  where  the  junction’s  small  capacitor-like  behavior  oscillates 
at  the  applied  frequency  [1],  [32J.  When  light  is  coupled  across  the  barrier 
causing  Fermi-level  modulation  at  this  optical  frequency,  the  second-order 
nonlinearities  of  the  l-V  curve  result  in  a  rectification  of  the  applied  optical 
voltage,  causing  a  DC  voltage  to  develop  across  the  barrier.  Since  this 
second-order  effect  is  proportional  to  the  second-order  derivative  of  the  l-V 
curve,  the  rectified  voltage  observed  versus  a  bias  voltage  will  trace  the 
second  derivative  of  the  l-V  curve.  These  second  derivative  features  appear 
as  resonant  curves  with  zero’s  at  the  inflection  points  of  the  l-V  curve  due 
to  the  band  gap  dictated  by  the  Fermi-Dirac  distribution  at  the  operating 
temperature.  In  other  words,  consider  impressing  a  sinusoidal  voltage 

F(t)  =  VpSin(Wj,t)  (2.4.10) 
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upon  the  DC  bias  voltage  being  applied  to  the  junction.  Assuming  that  the 
modulating  potential,  v^,  is  small  enough  to  produce  negligible  changes  in 
K^b)'  current  can  be  expressed  as  the  expansion 


(2.4.11) 


From  this  expression  it  is  obvious  that  if  the  detecting  apparatus  averages 
the  AC  tunneling  current,  then  only  the  even-powered  derivatives  will 
manifest  a  response  since  =  0,  while  <V(2")(t)>t  =  (1/2)Vo^.  Asa 

direct  consequence  of  this,  the  rectified  signal  will  appear  in  the  second 
derivative  term,  if  reotification  occurs,  yielding 


i 

4 


(2.4.12) 


Therefore,  the  second  derivative  of  the  junction’s  l-V  curve  can  be  plotted 
versus  V^,,  by  identifying  this  DC  rectified  response. 

Now,  assuming  a  barrier  shape  not  of  a  simple  square  barrier,  but  one  that 
embodies  an  asymmetry  such  as  that  in  Figure  7,  it  is  possible  to  determine 
,  and  L.  Several  references  [141[33]  describe  the  calculations  and  derivation 
leading  to  the  final  result  of 


[  ^ 

Sexp[S]  Sexp[S 

Sjcm^e 

L  "  j 

324q)a 

(2.4.13) 
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where  S=  1.025  L  \^,  a  is  the  junction  area  in  -  the  junction's 

impedance  --  is  in  Ohms,  L  (the  barrier  thickness)  is  in  Angstroms  and  9 
(the  average  barrier  height)  is  in  units  of  eV  (see  Figure  7).  From  this 
equation,  and  impedance  and  capacitance  measurements  (assuming  a  parallel 
plate  capacitor)  of  the  junction,  it  is  then  possible  to  find  9],  L,  and  an 
asymmetry  factor,  a,  where 


-  A(p  _  cp2  -  (Pj 
2^  2^ 


(2.4.14) 


once  <P2  has  been  determined  from  Fowler  plots.  Finally,  it  is  extremely 
important  to  remember  that  the  above  results  apply  only  when  the  impressed 
AC  signal  modulates  the  Fermi  level  of  one  electrode  --  i.e.,  the  AC  field  is 
coupled  across  the  barrier. 

2.4.3  Differential  Heating 

Finally,  of  all  the  current  producing  responses  the  remaining  major 
contributing  one  comes  about  as  a  result  of  the  way  the  experiment  was 
conducted.  As  can  be  seen  in  Figures  8  and  9,  by  using  a  phase  sensitive 
detection  technique  (thus  requiring  that  the  impinging  laser  light  on  the 
junctions  be  chopped)  resulted  in  a  slight  distortion  of  the  l-V  curves  due  to 
thermal  heating. 

Since  the  phase  sensitive  detection  technique  determines  a  signal’s 
presence  by  beating  a  reference  signal  against  the  incoming  signals  to  be 
analyzed,  thus  differentiating  the  intended  information  from  amidst  the 
noise,  it  is  possible  to  visualize  an  induced  differential  heating  signal 
sneaking  its  way  through  into  our  data.  The  unfortunate  part  about  this 
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signal,  at  first  inspection,  is  that  it  could  possibly  have  a  very  similar 
current  response  as  the  sought  after  rectified  signal. 

In  essence,  then,  what  the  lock-in  would  be  seeing  and  manipulating  is 

I(V,T+AT)-1(V,T)  (2.4.15) 

which  is  merely 

at  (2.4.16) 

dT 

That  is,  that  the  first  derivative  of  the  l-V  curve  with  respect  to 
temperature  could  possibly  contribute  to  the  rectified  signal.  Fortunately, 
as  will  be  shown  later,  once  the  temperature  dependence  of  the 
superconducting  band  gap  is  taken  into  account,  there  is  no  question  that  the 
observed  signal  is  a  result  of  rectification. 
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Differential  Heating 
Comparison  of  i-V  Curves 


CM 
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Figure  8  --  Illustration  of  the  Heating  Effect  on  the  Junction 
Due  to  Illumination 
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Figure  9  --  Calculated  Differential  Heating: 
[l(V)]illuminated  -  [l(V)]unilluminated 
(1mW,  5145A  ) 
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Chapter  3 

Experimental  Arrangement 

■  1  Junction  Fabrication^ 

Crucial  to  the  entire  project  was  the  quality  and  fabrication  of  the 
metal-oxide-metal/superconductor  tunneling  elements.  These  junctions 
were  made  via  a  metal  evaporation  process.  Basically,  the  quartz  substrate 
was  masked  by  a  thin  metal  plate  with  two  identical  slits  in  it,  diagonally 
offset  (see  Figure  10).  The  general  dimensions  of  these  slits  were  40pm 
wide  by  2.75mm  long.  It  was  important  to  ensure  that  these  masks  were  In 
intimate  contact  with  the  substrate,  so  as  to  avoid  "bleeding"  of  the 
evaporated  metal  underneath  the  mask. 


Figure  10  --  Substrate  Masks  for  Thin  Film  Deposition 


^For  an  in-depth  discussion  of  the  fomnation  of  metai-oxides  films,  see  [34],  [35]. 
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\ 

Once  the  masks  were  in  place,  the  junctions  were  formed  by  first 
evaporating  Al  to  a  height  of  100-1 50A,  measured  via  a  quartz  resonator.  (In 
addition,  during  the  earlier  attempts  to  refine  the  technique  of  depositing 
these  thin  films,  several  Sloan  Dektak  II  stylus  measurements  were  made  to 
confirm  lead  heights).  This  lead  was  then  oxidized  using  a  glow  discharge  of 
15-20  seconds  with  the  overall  conditions  centering  around  50p  of  O2.  while 
the  discharge  current  was  320  mA.  This  typically  produced  an  oxide 
thickness  of  about  10-20A.  Without  exposing  the  freshly  oxidized  surface  to 
the  atmosphere,  the  masks  were  rotated,  and  the  Vanadium-Titanium  alloy 
(50%-50%)  was  then  deposited  (again  by  evaporation)  at  right  angles  to  the 
two  Al  strips,  to  a  height  of  500A  (Figure  11).  Once  the  VTi  strip  had  been 
deposited,  a  new  mask  with  a  square  hole  was  then  placed  on  top  of  the 
junction,  allowing  a  protective  coating  of  500A  of  MgO  to  be  deposited.  This 
protected  the  critical  area  of  the  junction  from  further  oxidation,  moisture, 
and  physical  erosion.  This  completed  the  basic  structure  of  the  tunneling 
elements. 


Figure  11  -  The  Resulting  AI-AI  oxide-VTi  Tunneling  Element 
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As  for  the  "care  and  feeding"  of  the  junctions,  any  long  term  preservation 
was  accomplished  by  keeping  them  in  a  desicator,  pumped  down  as  far  as  the 
rough  pumps  would  take  it,  which  was  approximately  30-40p  If  these 
elements  were  exposed  for  too  along  a  time  to  high  humidity  or  atmosphere, 
there  was  a  tendency  for  them  to  increase  in  impedance,  sometimes 
reversible  by  vacuum  pumping  the  samples. 

For  historical  reasons,  if  for  no  other,  it  should  be  mentioned  that  for  a 
significant  amount  of  time,  we  attempted  to  use  AI-AI  oxide-Pb  as  the 
Junction  materials.  In  addition,  the  initial  attempts  were  with  3  pm  wide 
leads.  At  these  dimensions,  and  with  the  intended  film  thicknesses,  it 
appears  as  if  Pb's  amorphous  behavior  was  the  cause  of  the  many  difficulties 
experienced.  Not  only  does  lead  tend  to  migrate,  but  at  these  dimensions, 
dendrite  formations  tended  to  "puncture"  the  oxide  layer,  resulting  in  very 
rapid  deterioration  of  the  junction  impedances.  Junction  lifetimes  were 
between  one  and  three  days  at  room  temperature,  and  often  the  thermal 
contraction  that  resulted  by  cooling  to  liquid  Helium  temperatures 
accelerated  this  "puncturing"  deterioration.  As  a  result,  we  sought  a  more 
stable  superconducting  metal,  and  found  VTi.  Possible  adaptive  solutions  are 
noted  in  Chapter  5. 

In  order  to  electrically  connect  the  junctions  to  the  diagnostics,  36  gauge 
copper  magnet  wire  was  silver  painted  to  the  large  square  pads  at  the  ends 
of  the  metal  leads.  The  other  ends  of  these  magnet  wires  were  then  hard 
soldered  into  a  printed  circuit  card  that  had  a  hole  cut  in  the  middle  of  it  (to 
allow  the  laser  light  to  strike  the  target  area).  This  printed  circuit  card 
then  could  be  inserted  into  a  connector  at  the  base  of  the  sample  rod  in  the 
cryogenic  dewar. 
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During  all  room  temperature  work,  great  attention  was  paid  to  the 
current  densities  that  the  leads  may  be  subjected  to.  Extreme  care  was 
exercised  to  avoid  exposing  especially  the  Aluminum  lead  to  anything  greater 
than  1000  Amps/cm2  . 

Finally,  in  comparing  room  temperature  impedance  measurements  with 
those  impedance  measurements  done  at  liquid  helium  temperatures,  normally 
the  junction  exhibited  no  more  than  a  10-20%  impedance  change.  Anything 
more  than  that  typically  indicated  a  poor  silver  paint  attachment. 

3.2  Cryogenic  Apparatus  and  Optics 

Most  of  the  experiments  were  performed  while  the  junctions  were  at  or 
below  liquid  Helium's  boiling  point  of  4.2“K.  All  the  experiments  done  at 
cryogenic  temperatures  involving  laser  irradiation  were  performed  at  or 
below  the  lambda  point  of  2.18'‘K  to  ensure  optical  clarity.  The  dewar  used 
in  this  research  was  a  product  of  Andonian,  having  a  liquid  Helium/sample 
chamber  volume  of  3  liters,  surrounded  by  a  vacuum  space,  which  is  again 
surrounded  by  a  liquid  Nitrogen  thermal  shroud  chamber,  which  is  finally 
surrounded  by  yet  another  vacuum  space.  All  vacuum  spaces  were  evacuated 
to  about  10’^  Torr  using  a  diffusion  pump.  The  Andonian  dewar  was  fitted 
with  an  extension  that  tailed  about  12  inches  below  the  3  liter  section.  It  is 
Into  this  tail  section  that  the  sample  rod  with  the  substrate,  junction  and 
printed  circuit  card  were  inserted.  Once  the  liquid  Helium  was  added,  a  total 
experimental  run  time  of  about  6-7  hours  could  be  expected,  prior  to  "burn¬ 
out"  of  the  liquid  helium  (this  is  to  include  pumping  on  the  He  section  with  a 
rough  pump  to  reduce  the  temperature).  All  temperature  measurements  were 
made  using  a  Lake  Shore  silicon  diode. 


35 


Once  in  place,  the  junction  could  be  targeted  with  the  lasers  through 
either  of  two  pairs  of  quartz  windows  (two  windows  in  the  front,  and  two  in 
the  rear),  each  pair  consisting  of  one  room  temp  "warm"  window  between 
atmosphere  and  the  vacuum  chamber,  and  a  liquid  helium  "cold"  window, 
between  the  vacuum  chamber  and  the  sample  area.  The  sample  rod,  and  thus 
the  substrate,  could  be  rotated  through  360“;  and,  with  a  micrometer  affixed 
to  the  top  of  the  sample  rod,  the  height  (relative  to  the  optical  plane)  could 
be  adjusted.  Horizontal  alignment  and  adjustment  was  obtained  by  moving 
the  final  mirror  and  final  focusing  lens  in  unison  (Figures  12  and  13). 

Using  an  Argon-ion  laser  that  was  displaced  from  our  lab  necessitated 
the  use  of  a  fiber  optic  cable  In  order  that  we  may  have  access  to  the  needed 
light.  The  wavelengths  used  were  the  six  strong  lines  between  4579A  and 
5145A.  The  beam  profile  was  gaussian,  ensuring  a  minimum  focal  spot  on  the 
junction  —  approximately  3-5  pm.  To  avoid  intra-lens  interference,  all 
lenses  were  AR  coated.  Focusing  of  the  beam  onto  the  junction  and 
"navigating"  the  focal  spot  around  junction’s  area  was  accomplished  using  a 
combination  of  techniques.  Several  direct  optical  viewing  options  were  tried 
to  no  great  level  of  success.  Finally,  by  simple  use  of  a  clean,  white  "screen" 
placed  behind  the  tail  section,  one  could  view  the  transmitted  light  and  the 
resulting  junction's  shadow.  It  was  easy  to  determine  from  this  shadow's 
image  when  the  focal  spot  was  not  on  the  junction,  and  approximately 
whether  or  not  the  beam  waist  was  actually  on  the  junction.  At  this  point,  all 
adjustments  were  performed  with  reference  to  the  electrical  signals 
obtained. 

Finally,  a  pair  of  polarizers  were  placed  in  the  beam  path  (after  the  pin 
diode  sensing  point)  in  order  that  the  power  and  polarization  could  be 
adjusted.  All  power  readings  were  taken  CW  after  the  polarizers  and  prior  to 
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the  first  of  the  two  final  mirrors,  and  the  correlation  factor  was  determined 
for  the  corresponding  power  at  the  junction's  location.  The  data  were 
corrected  using  this  calibration. 

3.3  Photoemissive  Barrier  Height  Determination  Arrangement ° 

Figure  12  outlines  the  experimental  arrangement  for  the  oxide's  energy 
barrier  height  determination.  These  data  were  obtained  strictly  at  room 
temperature,  and  as  a  result  the  thermal  contributions  had  to  be  taken  into 
account.  The  output  of  the  Ar"^  laser  was  tuned  to  a  single  line,  and  was  then 
focused  onto  the  junction  (again,  all  power  measurements  were  taken  after 
the  polarizers  and  while  CW).  A  mechanical  chopper  was  placed  in  the  beam 
path  at  the  focal  point  of  two  equal  achromatic  planoconvex  lenses,  providing 
a  70  psec  pulse  with  a  rise  time  of  3  psec,  and  an  illumination  rep  rate  of  38 
Hz.  The  junction's  response  at  the  beginning  of  the  pulse  was  determined  by 
extrapolating  the  signals  taken  at  15  psec  and  60  psec  after  the  pulse  was 
initiated  by  using  a  two-channel  boxcar  with  a  1  psec  gate  width  set  to 
average  300  samples.  The  pulse  initiation  was  determined  by  a  Thor  Labs 
silicon  pin  diode  which  was  impedance  matched  to  ensure  a  faithful 
representation  of  the  optical  pulse  shape.  This  arrangement  allowed  the 
determination,  and  thus  the  subtraction,  of  the  slow  laser-induced  thermal 
response.  The  resulting  power  normalized  data  as  a  function  of  the  bias 
voltage  provided  the  information  for  a  Fowler  plot,  yielding  the  barrier 
height  (Figure  17). 

As  was  determined  by  Elchinger  [1],  the  examined  largest  changes  in  L  and 
9  that  could  occur  through  the  drastic  temperature  change  from  300“K  to 
liquid  Helium  are  smali  (in  comparison  with  measurement  errors).  Reasonabiy 

^For  the  most  part,  this  procedure  faithfully  follows  that  which  is  described  in  [32]. 
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small  changes  in  R^j  (no  more  than  20%)  with  temperature  were  ignored  when 
determining  the  barrier  parameters. 

3.4  l-V  and  l"-V  Arrangement 

Figure  13  outlines  the  experimental  arrangement  for  the  l-V  and  l"-V  data 
collection.  All  data  were  collected  below  2.18*K,  with  the  stronger  signals 
being  obtained  at  I.T'K.  Again,  the  Ar"*^  light  was  tuned  to  a  single  line  and 
focused  onto  the  junction.  This  time,  though,  a  chopper  blade  was  used  that 
chopped  the  beam  at  400  Hz,  providing  a  pulse  of  1.25  msec,  while  the  pin 
diode  signal  was  used  as  the  lock-in  reference. 

For  both  the  barrier  height  determinations  as  well  as  the  l-V,  and  r'-V 
work,  the  electronics  design  was  basically  the  same,  requiring  a  bias  supply, 
a  lock-in  amplifier,  and  a  display.  For  the  barrier  height  work,  there  were  an 
additional  two  stages  of  amplification  that  the  l-V  work  did  not  have  (see 
Figure  12).  In  short,  the  junction's  response  was  first  amplified  by  a 
Tektronics  1A7A  differential  amplifier,  which  was  then  fed  to  the  boxcar. 
Aside  from  that,  the  two  arrangements  were  basically  the  same. 

The  0-1 50V  bias  voltage  supply  was  first  fed  into  the  Bias  Box  (Figure 
14),  which  incorporated  an  operational  amplifier  driven  feedback  loop  to 
ensure  voltage  control  when  the  junction  resistance  changed  dramatically 
with  bias.  In  addition,  a  series  current  limiting  resistor  was  included  at  the 
output.  This  resistor  was  made  variable  so  as  to  obtain  the  maximum  signal 
for  a  given  bias  voitage.  Reference  36  develops  the  second  derivative 
detection  circuitry  that  was  basically  used  in  this  experiment,  and  it 
indicates  that  the  ideal  circuit  would  have  the  bias  resistance  very  large  in 
comparison  with  R^,  effectively  making  a  current  supply  of  the  Bias  Box. 
However,  the  resulting  bias  developed  across  the  junction  is  very  small.  At 


39 


the  other  extreme,  if  the  bias  resistance  is  too  small,  it  ends  up  loading  the 
junction,  again  producing  a  small  response.  Consequently,  the  variable 
arrangement  in  Figure  14  was  used.  In  both  arrangements,  a  dual  staged 
amplifier  (Figure  15)  was  used  to  increase  the  overall  signal  strength  and  to 
minimize  the  digitization  error  prior  to  entry  into  the  lock-in. 

So  as  to  ensure  accurate  measurements  of  the  junction  impedance  both  at 
room  temperature  and  at  liquid  Helium  temperature,  4-pole  measurements 
were  made  (where  possible).  This  permitted,  in  essence,  cancellation  of  the 
inherent  line  impedance,  so  that  in  effect  only  the  junction  was  being 
examined  (see  appendix  3). 

Finally,  electrical  noise  was  a  constant  concern.  An  extreme  amount  of 
time  and  effort  was  expended  to  minimize  noise.  Great  care  was  taken  to 
avoid  ground  loops,  establishing  a  common  instrumental  ground  at  the  Bias 
Box.  All  wires  were  shielded,  including  the  junction  (by  the  dewar).  Low  band 
pass  filters  were  used  in  the  Amplifier  Box  (Figure  14)  to  help  filter  any 
unwanted  signals  prior  to  amplification.  And,  finally,  to  help  filter  out  some 
extremely  high  freq  noise  that  we  just  could  not  get  rid  of,  some  low  pass  n- 
filters  were  used  at  the  lock-in’s  input  which  significantly  cleaned  our 
signal. 
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Figure  13  --  Optical  Rectification  Experimental  Set-Up 
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Chapter  4 

Metal-Barrier-Superconductor  Response 

4.1  Barrier  Height  Determination 

Several  investigators  [1][29][31][42]  have  extensively  studied  tunneling 
barrier  shapes,  having  begun  with  trapezoidal  assun^ptions.  In  these  works 
the  barrier  height  (see  Figure  7)  of  %  has  been  measured  by  photoemission 
over  the  barrier  using  the  zero  current  intersection  in  a  Fowler  plot  as  the 
value  of  that  height.  In  addition,  most  often  the  q)j  barrier  height  could  also 
be  determined  if  a  large  enough  positive  bias  was  applied  to  the  metal  on 
that  side.  Due  to  the  dimensions  used  in  this  work,  the  necessary  bias  would 
have  resulted  in  a  current  density  in  the  Aluminum  leads  that  would  have 
been  risky  at  best,  and  damaging  at  worst.  Consequently,  this  necessitated 
calculating  the  parameters  of  L  and  (p,  from  capacitance  measurements,  once 
(Pj  was  known,  using  equation  (2.4.13),  and  the  definition  ^  ^  ^2 

By  observing  the  photo-induced  tunneling  currents  (PITC),  cpj  was 
determined  to  be  1.86  eV.  Figure  16  shows  a  typical  junction's  room 
temperature  PITC  response  to  several  wavelengths  of  light.  From  this  data, 
it  is  then  possible  to  take  VJpg  and  plot  it  against  hv.  The  resulting  Fowler 
plot  of  Figure  17  permits  a  best  straight  line  fit  extrapolating  to  the  Jp^  =0 
point  to  determine  (p2. 

Assuming  that  the  junction  forms  a  parallel  plate  capacitor,  then  from  a 
capacitance  measurement  and 

R^Cd=!!^R^  =  8.85x  (4.1.1) 
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Figure  1 7  --  Fowler  Plot  to  Determine  Barrier  Height 
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it  is  possible  to  determine  L.  Here,  these  units  permit  R^j  to  be  in  ohms,  L  to 
be  in  A,  and  a  to  be  in  pm2.  The  dielectric  constant,  £,  has  been  studied  by 
several  investigators  and  has  been  determined  to  be  8  [1][18][31][35][42].  For 
resistances  on  the  order  of  about  7  kohm,  the  capacitance  was  between  70- 
80  pF.  This  typically  implied  an  oxide  layer  thickness  of  about  14.5  A. 

Through  equation  (2.4.13),  and  with  L  and  cpj  determined,  this  leads  to  ^  = 
1.74  eV  and  thus  %  =  1.62  For  the  most  part,  then,  the  trapezoidal  tunneling 
barrier  has  been  characterized. 

4.2  Photoemission  and  Thermal  Responses 

All  of  the  above  barrier  parameter  determinations  were  accomplished  at 
room  temperature.  From  Figure  16,  it  is  possible  to  see  that  there  is  a 
sizable  PITC.  In  fact,  at  this  ambient  temperature,  it  is  the  PITC  that 
dominates  as  can  be  seen  by  observing  Figure  18,  the  70  psec  laser  pulse  at  a 
bias  voltage  of  -250  mV.  The  sloping  top  of  this  pulse  is  due  to  the  thermal 
heating  of  the  junction.  Obviously  here,  the  thermal  contribution  to  the 
entire  signal  is  easily  dominated  by  the  PITC. 

However,  once  these  junctions  are  cooled  to  liquid  Helium  temperatures, 
the  PITC  becomes  less  significant,  and  is  easily  dominated  by  the  thermal 
signal.  This  is  quite  evident  in  Figure  19:  by  illuminating  the  junction  with 
1  mW  of  5145A  light,  it  is  possible  to  see  that  at  the  higher  bias  voltages, 
the  signal  is  very  nearly  zero  (thus  very  little  PITC  contribution).  In  addition, 
since  there  is  no  point  at  which  the  curve  intersects  the  x-axis  (except  for 
the  origin),  there  are  no  inflection  points  indicated  and,  thus,  at  this  power 
setting,  the  thermal  contributions  also  dominate  the  rectified  response. 

Not  only  does  this  curve  nicely  match  the  "calculated"  thermal  response 
curve  of  Figure  9,  but  the  theoretically  predicted  thermal  response  curves  of 
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Figure  18  --  Fast  Photoresponse  with  Thermal 
Response  Superimposed  at  -250  mV  Bias 
(T  =  300'K) 
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Curves  demonstrating  the  effect 
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Calculated  ‘Differential  Thermal  signals 
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Figure  22  --  dl(V)/(lT  of  theoretical  l-V  curves, 

Showing  that  the  infinitesmal  differential  temperature  contributions 
show  no  inflection  points  in  l-V. 
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Figure  20  and  Figure  21  support  the  general  trace  shape  as  well.  The  signal 

arising  from  the  change  in  temperature  of  the  metal-barrier-superconductor 

caused  by  laser  illumination  is  clearly  different  from  the  predicted  rectified 

response  (Figure  30).  Giaever  [17]  showed  that  upon  heating  a 

superconducting  junction,  the  l-V  curve  smoothly  changes  toward  a  straight 

line,  the  slope  of  which  represents  the  normal  junction's  high  bias 

conductivity.  Since  the  higher  temperature  curve{s)  (Figure  20)  never 

intersect  the  unilluminated  ,  iower  temperature  curve,  there  can  be  no 

response  that  changes  sign  (Figure  21  and  Figure  22).  What  this  indicates, 

then,  is  that  the  thermal  response  signal  can  never  demonstrate  that  the 

original  l-V  curve  has  an  inflection  point  and,  thus,  is  fundamentally 

different  from  the  rectified  response.  In  addition,  this  does  imply,  though, 

that  the  thermal  contributions  will  alter  the  relative  peak  ratios  from  that 

predicted  by  theory.  To  put  a  final  nail  in  the  thermal  coffin,  all  of  the  above 

thermal  examinations  have  been  carried  out  on  finite  temperature  changes. 

Figure  22  displays  the  ^ curve  (for  an  infinitesimal  temperature  change) 

dT 

of  equation  (2.3.5)  taking  into  account  the  temperature  dependence  of  the 
superconducting  band  gap  via  the  BCS  relation.  Clearly,  this  shows  no 
inflection  point  as  well. 

However,  this  does  not  yet  explain  why  the  photoemissive  signal  has 
vanished  at  the  liquid  helium  temperatures.  A  well  known  characteristic  of 
metals  is  that  when  the  metal  is  cooled,  its  reflectivity  increases.  This  is 
due  to  a  reduction  in  the  effective  skin  depth.  For  good  conductors,  the  skin 
depth,  6,  can  be  written  [43] 


6=  g  - 
42n  p  fo  o 


(4.2.1) 
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where  c  is  the  speed  of  light,  p  Is  the  permeability,  to  is  the  circular 
frequency  of  the  impinging  light,  and  o  is  the  conductivity  (Gaussian  units). 
From  the  data  gathered  in  references  [44-46]  on  the  room  and  low 
temperature  electrical  conductivity  of  Aluminum,  equation  4.2.1  yields 


^300»K 

^2‘IC 


=  64A 


=  6  A 


(4.2.2) 


Assuming  that  ^«1,  it  is  possible  to  show  that  the  coefficient  of 

A. 

"penetrating"  light  (that  light  which  is  not  reflected,  but  in  some  manner 
scatters  or  interacts  with  the  Aluminum)  is 


fans 


(l-R) 


6  4jtn, 


(4.2.3) 


where  n-i  is  the  index  of  refraction  for  the  medium  from  which  the  light 
penetrates  the  metal  (in  this  case,  either  superfluid  Helium,  or  a  vacuum). 

Finally,  the  mean  free  path  length  of  the  photoelectrons  is  vastly 
different  at  the  two  temperatures.  This  factor,  too,  will  effect  the  overall 
amplitude  of  the  photoemissive  current.  Chaverri,  et.  al.,  [47]  have  studied 
Aluminum  both  at  room  temperature  as  well  as  at  low  temperatures.  From 
their  data,  and  using  the  Drude  model  to  extrapolate  to  2°K,  the 
photoelectron  mean  free  path  lengths  at  the  respective  temperatures  are 

,  (4.2.4) 

Ar2„K  =950  A 
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Now,  from  equation  (2.4.9),  and  from  reference  [29],  the  photoemissive 
current  can  be  written  as 

■pE-  22^(hv-(P2f(ijf  (4.2.5) 

where  P  is  the  power  of  the  impinging  light.  Expanding  f, 

■p.E- 2^(hv-<P2f(j^)o^(i-R)  H-Z-e) 

where  o  is  the  cross  section  for  photogeneration  of  hot  electrons,  and  x  is 
the  mean  free  time  of  photoexcited  electrons.  Rewriting  the  above  equation 
in  terms  of  mean  free  path  lengths,  and  Fermi  velocities,  Vf, 

where  a  is  the  light  absorption  coefficient,  and  is  the  number  of  electrons 
Illuminated.  Finally,  this  can  be  rewritten  in  terms  of  the  skin  depth,  6p  and 
the  electron  density,  6„  ,  as 

IpB.-  2fP(hv-<p,fW(a)!5  yjvolume  illuminated)  |6„  Ji  |^|  (1-R)  (4.2.8) 

'tlV  /  \  f  f 

Measurements  were  made  at  both  room  temperature  and  at  superfluid  Helium 
temperatures  at  the  extreme  ends  of  the  bias  voltage  ramp.  It  was  found, 
using  the  same  light  intensity  at  both  temperatures,  that  the  overall 
reduction  in  photoemissive  signal  strength  was  reduced  by  a  factor  of  about 
20  when  at  liquid  helium  temperatures. 
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Now  assuming  that  a,  VpCpj  and  the  probability  for  tunneling  remain 


unchanged  between  the  two  temperatures,  than  the  relative  photocurrent 
intensities  will  be 


^°K  ^  ( 1  R2°k)  ^2°»C  ^2°K. 

^300°K  (  ^  ~  ^300‘‘K  ^300"K 


JL 

10 


(4.2.9) 


which  is  not  a  bad  correlation  with  the  overall  reduction  seen  here.  And, 
considering  that  the  photoemissive  contributions  are  more  intense  at  the 
higher  bias  voltages,  this  bodes  well  for  the  lack  of  any  visible  contribution 
from  the  photoemissive  current  during  the  rectification  studies,  since  these 
data  were  collected  at  voltages  near  zero  bias  {±5  mV). 

4.3  Fermi  Level  Modulation  at  Optical  Frequencies  --  Rectification  Data  and 
Response 

Armed  with  the  above  appreciation,  it  is  now  possible  to  show  that  when 

illuminating  the  metal-barrier-superconductor  junction  with  the  proper  light 

intensity  ,  the  responses  seen  are  dominated  by  rectification  due  to  the 

modulation  of  the  Fermi  level  at  an  optical  frequency. 

2 

From  equation  (2.3.5),  — has  been  calculated  and  displayed  in  Figure 

dV^ 

23.  Again,  note  here  that  the  significant  qualitative  difference  between  the 
rectified  and  thermal  signals  is  that  the  thermal  response  does  not  change 
sign  and  thus  does  not  indicate  an  inflection  in  the  l-V  curve,  in  comparison 
to  the  rectified  response  which  does.  Both  of  the  curves  in  Figure  23  typify 
the  metal-barrier-superconductor  (M-B-S)  junction  response  because  of  the 
familiar  quick  rise  in  current  at  Vj,  =  A/e,  and  the  smooth,  asymptotic 
behavior  toward  linearity  at  large  bias.  From  this  behavior,  it  is  also  to  be 
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Figure  24  -  T-V  of  a  S-B-S  Junction 
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(Kudio  Second  Derivative 
modulation,  2mV  peak*to-peak) 


determined  that  the  junctions’  responses  were  truiy  of  a  normai  metai- 
barrier-superconductor  response.  As  Giaever  [17]  so  aptly  demonstrates,  the 
l-V  and,  thus  the  l"{V)  curves  are  vastly  different  when  the  junction  is 
comprised  of  a  superconductor-barrier-superconductor.  Careful  monitoring 
of  the  second  derivative  response  combined  with  ensuring  that  the 
temperature  was  no  colder  than  1.7*K  guaranteed  that  Aluminum  remained 
normal.  In  addition,  impedance  measurements  on  Aluminum  at  the  operating 
temperatures  always  indicated  a  low  but  finite  (approximately  5-15  Cl) 
resistance.  Finally,  to  illustrate  the  fundamental  difference  between  a  M-B- 
S  and  an  S-B-S  rectified  response.  Figure  24  clearly  displays  a  zero  slope 
region  inside  2A  —  consistent  with  the  second  derivative  of  an  l-V  curve  for 
two  superconductors  who  have  no  current  flowing  while  the  respective  band 
gaps  are  "juxtaposed".  Consequently,  it  was  easy  to  ensure  that  the  junction 
was  truly  M-B-S  during  the  experiments. 

Now  in  section  4.2,  it  was  shown  that  although  the  thermal  response  Is 
dominated  by  the  rectified  response  (that  is,  the  rectified  response  changes 
sign  while  the  thermal  response  does  not),  the  thermal  contributions  can 
alter  the  relative  peak  ratios.  A  ready  means  of  demonstrating  this  is  by 
applying  an  audio  frequency  across  the  junction.  Figure  25  shows  the 
rectified  response  of  a  M-B-S  junction  when  a  2  mV  (peak-to-peak),  200  Hz 
sinusoidal  signal  was  applied.  The  resulting  heating  effect  is  minimal  (but 
still  observable),  while  the  trace  is  very  nearly  that  predicted  by  theory. 

Note  here,  also,  that  the  height  of  the  peaks  (A  and  B)  should  be  affected 
by  the  temperature,  through  the  temperature  dependence  of  the  l-V  curve,  as 
well  as  through  V^gipevyave-  Since  the  rectified  response  is  linear  in 
V^sinewave-  the  ratio  of  A/B  is  not  a  function  of  Vsinewave-  t)ut  primarily  of 
temperature. 
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This  can  be  seen  in  the  series  of  Figures  26-29,  as  the  power  of  the  iaser 
light  was  changed  from  120  pW  to  20  pW.  The  ratio  of  A/B  remains  nearly 
constant,  approximately  18:1,  in  this  power  region  where  the  heating  effect 
had  been  minimized.  In  addition,  Figure  30  shows  a  theoretical  superposition 

of  the  — curve  on  top  of  the  ^  curve.  Note  here  that  the  relative 
dT 

peak  heights  approach  a  ratio  of  A/B  =  10:1.  Clearly  this  indicates  that  the 
differential  thermal  signal  is  playing  a  part,  although  an  additional  factor  of 
2  (or  more)  needs  to  be  accounted  for. 

Elchinger  [1]  accomplished  an  exhaustive  study  of  the  response 
dependence  upon  polarization  and  the  angle  of  incidence.  While  aware  of  this 
dependence,  we  only  maximized  the  response  amplitude  by  ensuring  that 
there  was  an  E-field  perpendicular  to  the  junction.  However,  this  point  is 
worthy  of  some  discussion  for  it  has  been  the  grounds  for  discounting 
rectification  via  modulation  of  the  Fermi  level  at  optical  frequencies  [1]. 
Elchinger  notes  that  the  absorption  is  Independent  of  polarization  at  normal 
Incidence,  while  for  increasing  incident  angle  the  absorption  decreases  when 
the  electric  vector  is  perpendicular  to  the  pane  of  incidence.  When  the 
polarization  is  parallel  to  the  incident  plane  the  absorption  slightly 
increases  until  the  incident  angle  reaches  the  principle  angle.  Although 
Elchinger  observed  the  characteristic  response  at  normal  incidence  (and  with 
a  large  signal  response),  he  attributes  the  fact  that  since  there  is  no  E-field 
component  perpendicular  to  the  barrier  that  Fermi  level  modulation  cannot 
occur.  To  this  end,  on  an  optical  scale,  the  surface  of  the  target  area  cannot 
be  perfectly  smooth  and  thus  cannot  be  completely  normal  to  the  incoming 
light.  Hence,  there  must  be  some  component  of  the  E-field  which  will  be 
perpendicular  to  the  barrier. 
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Finally,  some  discussion  is  warranted  regarding  the  overall  amplitude  of 
the  rectified  response  that  is  observed.  Taking  Figure  27  as  an  example, 
focusing  the  light  to  an  area  of  5  pm^,  the  resulting  square  of  the  amplitude 
of  the  E-field  would  be 


£2=87LP  =  208(^f 

ac  \cm/ 

Using  this  value,  equation  (2.4.12),  and  the  theoretical  value  for 


(4.3.1) 


d'  1(V) 

dV^ 


from 


Figure  23  results  in  an  expected  26  mAmps  across  the  7.2  kQ  junction. 

Hardly  the  signal  observed.  However,  taking  into  account  the  skin  depth,  the 
lead  thickness  (between  lOOA  and  150A),  and  the  resulting  attenuation  of  the 
E-field  that  goes  as  e“*^,  then  factoring  in  the  dielectric  constant,  a  resulting 
current  on  the  order  of  0.4  nAmps  could  be  expected.  This  is  very  nicely  in 
the  region  that  we  were  observing  signals,  indicating  that  the  E-fieid  was 
penetrating  the  Aluminum,  albeit  inefficiently. 

In  addition,  it  should  be  noted  that  using  the  conditions  reported  by 
Elchinger,  it  is  possible  to  very  nearly  exactly  theoretically  model  his 
rectified  response.  Assuming  no  heating  effect  at  all,  his  Ai-Al205-Pb 
M-B-S  junction  would  have  yielded  a  peak  ratio  of  3.7:1  (Figure  31). 

However,  he  reports  a  peak  ratio  of  7.6:1,  which  is  nicely  modeled  by 
accounting  for  a  heating  effect  due  to  illumination,  and  taking  into 
consideration  the  band  gap  dependence  on  temperature  (Figure  32).  The 
significant  difference  in  Elchinger's  conditions  and  those  in  this  work  is  that 
his  Aluminum  leads  were  4  to  6  times  thicker.  There  are  indications  [50]  that 
as  the  Aluminum  thickness  diminishes  to  dimensions  close  to  ours  that  the 
E-field  is  able  to  propagate  longitudinally  through  the  electron  plasma.  This 
may  have  the  resulting  effect  of  skewing  the  peak  ratios  as  we  have  seen. 
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Chapter  5 

Future  Work 

As  was  demonstrated  in  the  last  section,  rectification  as  a  result  of 
Fermi  level  modulation  across  the  barrier  dominates  when  the  E-field 
couples  across  the  barrier  -  resulting  in  a  response  that  traces  the  second 
derivative  of  the  junction’s  l-V  curve.  However,  the  strength  of  this 
rectification  signal  in  this  work  was  several  orders  of  magnitude  weaker 
than  desired  or  initially  expected  as  a  result  of  the  means  by  which  the  E- 
field  does  indeed  get  coupled  across  the  barrier.  That  is  to  say,  once  again, 
that  the  penetration  of  the  illumination  through  the  metal  lead  is  a  most 
inefficient  means  of  accomplishing  the  task  at  hand.  Better  signals  may  be 
obtained  through  some  clever  and  judicious  techniques. 

With  the  present  state-of-the-art  photolithographic  fabrication 
techniques,  it  is  possible  to  construct  high-speed  junctions  whose 
dimensions  are  <  1pm.  In  fact,  with  current  nano-architecture  capabilities, 
resolutions  down  to  300nm  are  possible.  At  these  dimensions,  it  is  possible 
to  see  from  equation  (4.1.1)  that  the  speed  of  this  element  will  be 
significantly  enhanced.  Using  any  of  a  variety  of  geometries,  the  aim  of  the 
fabrication  would  be  to  integrate  the  high-speed  junctions  with  optical 
frequency  antennae  resonating  as  a  dipole.  These  antennae  could  be  designed 
to  match  the  junction  impedance  (-100Q)  and,  as  a  result,  would  permit  the 
coupling  of  the  impinging  E-field  across  the  barrier  at  an  amplitude  sizably 
larger  than  the  E-field  produced  by  the  edge/penetration  coupling  technique 
used  in  this  work.  Not  only  would  this  permit  huge  signals,  and  great  signal- 
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to-noise  ratios,  but  would  also  allow  nanowatt  power  levels  of  light  to  be 
used.  This  is  not  only  a  distinct  improvement  over  the  resulting  current  seen 
in  this  work,  but  an  immense  improvement  over  the  earlier  works  depending 
on  the  relatively  minor  non-linearities  (about  five  orders  of  magnitude  less 
than  metal-oxide-superconductor)  of  metal-oxide-metal  elastic  tunneling 
elements  used  in  previous  works. 

At  this  point,  one  would  have  a  high-speed  electron  tunneling  element 
with  a  response  time  fast  enough  to  permit  me  flow  of  optical  current 
through  the  element  at  visible  and  UV  frequencies.  As  a  result  of  this,  and 
again  using  the  nonlinear  current-voltage  characteristics  of  the  metal- 
oxide-superconductor  tunneling  element  which  permits  it  to  be  a  broadband 
frequency  mixer/multiplier,  it  is  then  capable  of  being  used  to  directly  mix 
the  frequencies  of  lasers  differing  by  several  octaves  -  the  frequency  of  a 
microwave  standard  (or  klystron)  with  that  of  an  IR  laser,  the  frequency  of 
an  IR  laser  with  that  of  a  visible  laser,  or  even  a  UV  laser  -  thus 
eliminating  the  current  multi-stage  arrangement. 

With  the  broadband  frequency  mixing  capabilities,  this  extremely  high¬ 
speed  element  makes  it  possible  to  construct  a  phase-locked  frequency 
multiplier  chain,  forming  a  "transfer  oscillator".  It  is  this  transfer  oscillator 
that  would  consist  of  several  phase-locked  stages,  thus  extending  the 
microwave  frequency  measurement  technology  well  into  the  visible  and  even 
into  the  UV  region.  The  immediate  result  of  this  extension  will  make  it 
possible  to  measure  the  absolute  frequencies  of  the  visible  and  UV  lasers 
used. 

This  has  some  severe,  important,  and  beneficial  implications  for  the 
science  of  spectroscopy,  as  well  as  for  the  time  and  length  standards.  In 
addition,  it  will  now  be  possible  to  construct  a  master-clock  oscillator 
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operating  at  an  optical  frequency  with  an  output  in  the  microwave  region,  the 
microwave  output  being  phase-locked  to  the  master  optical  clock.  Or, 
conversely,  it  will  be  possible  to  construct  a  Cesium  clock  with  an  output  at 
an  optical  frequency,  with  this  optical  output  phase-locked  to  the  Cesium 
clock  frequency. 

Throughout  the  above  discussion  regarding  the  future  development  of  a 
true  master-oscillator,  and  hence  an  optical  clock,  it  must  be  realized  that 
what  is  required  is  a  laser  with  a  highly  reproducible  and  stable  frequency. 
What  makes  it  possible  to  have  such  a  stable  laser  is  the  ability  to  "time" 
optical  frequencies.  This  work  demonstrates  that  the  metal-oxide- 
superconductor  electron  tunneling  element  provides  this  abiiity  once  the 
Impinging  E-field  is  coupled  across  the  barrier,  resulting  in  the  Fermi  level 
being  modulated  at  the  applied  optical  frequency. 
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Appendices 

Appendix  1  -  Justification  of  Matrix  Eiement  Remaining 
Unchanged  when  the  Metai  Goes  Superconductingii 

An  expression  for  the  normal  to  superconducting  current,  l^s,  starting 
from  microscopic  theory  was  first  given  by  Cohen  et  al.  [37],  and  shortly 
thereafter  Schrieffer  [38-39]  took  a  slightly  modified  approach.  This 
explanation  follows  the  latter  article,  after  presenting  some  background 
information. 

In  1957,  Bardeen,  Schrieffer  and  Cooper  [40]  offered  the  first 
microscopic  theory  of  superconductivity.  The  main  assumptions  that  are 
critical  to  this  discussion  were; 

(a)  The  super  conducting  ground  state  can  be  expressed  solely  in  terms 
of  Cooper  pairs  so  that  the  states  (k,  -k)  are  occupied  or  empty 
simultaneously  (assuming  that  the  spins  are  and  remain  antiparallel  so  as  to 
have  no  need  for  spin  indices). 

(b)  The  various  interactions  may  be  taken  identical  in  the  normal  and 
superconducting  states  and  only  the  phonon  and  screened  Coulomb 
interactions  need  to  be  separated  for  attention. 

(c)  The  difference  between  the  phonon  and  screened  Coulomb 
interactions  is  -Vyj-  which  may  be  expressed  in  the  simple  form 

(V,  for  leJ,  IeJ  s  k0D 

Vh,  =  I  '  (A1.1) 

lO,  otherwise 

“I  '*  Background  information  taken  from  Ref[3]. 
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where  the  energy  is  measured  from  the  Fermi  surface  and  0d  is  the  Debye 
temperature.  Now,  admittediy,  this  is  an  obvious  oversimplification  which 
completely  disregards  the  details  of  the  interactions.  On  a  basic  level, 
though,  it  gives  a  good  approximation  because  the  superconducting  properties 
(in  reduced  coordinates)  are  hardly  dependent  on  crystal  structure  and  normal 
state  properties. 

Therefore,  in  accordance  with  the  above  assumptions,  the  ground  state 
energy  of  the  superconducting  state  at  T  =  0*K  (relative  to  the  energy  of  the 
normal  state)  may  be  written  as 

Wg  =  e,jV^  -  -2) 

where  is  the  probability  of  state  (k,  -k)  being  coupled  and  is  the 
probability  that  it  is  empty.  Obviously,  then, 

v^  +  u^=l  (A1.3) 


Now,  the  first  term  in  (A1.2)  gives  the  difference  in  kinetic  energy  between 
the  superconducting  and  normal  phases  at  T  =  O'K.  The  problem 
mathematically  is  to  minimize  the  ground  state  energy  as  expressed  in  (A1.2) 
with  respect  to  the  probability  v,j.  This  ends  up  being 


where 


(A1.4) 


(A1.5) 
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and  may  be  determined  from  an  integral  expression.'' 2 
Now,  Schrieffer  writes  the  matrix  element  in  the  form 

(A1.6) 

Thus,  the  matrix  element  is  not  precisely  identical  to  that  in  the  normal 
state,  since  there  is  this  additional  factor  of  .  In  addition,  another 
important  difference  is  that  there  are  two  quasiparticle  states  for  a  value  of 
energy  above  the  gap  as  shown  in  Figure  A1.  Thus  an  electron  of  momentum  kj 
on  the  left  side  may  tunnel  both  into  and  ^2- 


Figure  A1  -  Energy-Momentum  Diagram  of  an  M-B-S  Junction  at 
an  Applied  Voltage,  V,  in  the  Schrieffer  Quasiparticle  Representation 


“•2  Interestingly,  if  we  further  assume 

^  "  constant,  for  <  kBp 

^  |0,  for|EjJ>k9p 

and  A  «  kOp,  Nt4(£]()B  N(4(0)  (where  l^o)  is  the  density  of  states  at  the  Fermi  surface  in  the 
normal  state),  then  a  nicely  simple  expression  may  be  obtained  for  the  ground  state  energy 

Ws  =  -^Nn(o)a^ 
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Hence  the  current  will  be  proportional  to  the  sum 


w'K/>4)=rr|^2+^+y  (A1.7) 

But  according  to  {A1.5)  result  that  the  current  is 

proportional  to  |tJ^.  This,  then,  is  the  reason  for  the  correctness  of  the 

simple  tunneling  picture  and  the  phenomenological  theory  discussed  in 
chapter  2.3. 


Appendix  2  »  Faithful  Four  Pole  Impedance  Measurements 

The  technique  used  in  this  research  to  determine  the  impedance  of  the 
individual  leads,  as  well  as  the  impedance  of  the  junctions  was  via  "4-pole" 
measurements.  The  advantage  here  is  that  the  lead  resistances  of  the 
measuring  circuit  can  be  eliminated,  thus,  providing  faithful  readings  of  the 
object  being  measured. 

"4-pole"  measurements  are  conceptually  depicted  in  Figure  A2.1.  Across 
two  separate  leads,  a  current  (provided  by  the  Bias  Box  current  source)  is 
introduced  across  the  junction.  From  the  juxtaposing  quadrant,  voltage 
measurements  of  the  diode  are  then  taken.  It  is  this  geometry  of  the  circuit, 
combined  with  a  large  input  impedance  of  the  amplifiers  that  permit  the 
faithful  measurements. 


Figure  A2.1  -  "4-Pole"  Measurements  of  Junction  Impedances 

In  Figure  A2.2,  a  simplified  circuit  diagram  represents  the  experimental 
arrangement,  to  include  the  voltage  measurement  loop.  Rb  is  the  bias  resistor 
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in  the  Bias  Box.  The  various  Rl  identified  resistors  represent  the  inherent 
impedance  of  the  lines  going  to  and  coming  from  the  junction.  R^j  is  the 
representative  impedance  of  the  junction,  and  Ramp  represents  the  input 
impedance  of  the  opamp  in  the  Amplifier  Box. 


h 


Obviously  from  Kirchoff's  Current  Law  '^*0.  addition,  following 

a  standard  approach  for  a  current  divider,  it  is  possible  to  write 


i„  = 


^L1  ^L2 


°  (R^+Rl2)  +  (Ru  +  Ra  +  Ri2) 


(A2.1) 


Now,  rearranging  a  bit,  gives 


In  this  form  it  is  readily  apparent  what  happens,  since  the  fraction  in  the 
square  brackets  is  denominated  by  Ra  (which  is  on  the  order  of  108,  and 
typically  5  orders  of  magnitude  greater  than  R^).  Consequently, 


_  ^L2  j 

(Ru  +  Rl2+Ra) 


This  indicates  that  all  of  the  original  current,  ij,  flows  through  the  diode, 
thus  allowing  faithful  voltage  and,  therefore,  impedance  measurements 
across  it. 
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Appendix  3  »  Numerical  Integration  Using  the  MIT  Cray  X-MP 
Supercomputer 

In  order  to  calculate  the  l-V  and  T’-V  responses,  it  was  necessary  to 
perform  numerical  integrations  on  the  integral  in  equation  (2.3.5)  and  its 
associated  derivatives.  After  applying  two  coordinate  transformations, 
equation  2.3.5  loses  its  integratable  singularity,  and  is  placed  into  a 
computer  convenient  form.  Beginning  with  the  basic  form  of  the  above 
integral. 


I^3=AN„^0)N^0) 


E|  [f  (E-eV)  -  f  (E)l  dE 


E^-A^ 


1/2 


(A3.1) 


Integrating  for  all  possible  energies  yields 


*NS 


E[f(E-eV)-f(E)jdE 


E^-A^ 


1/2 


E[f(E-eV)-f(E)  dE 


E^-A^ 


1/2 


(A3.2) 


Introducing  x  +  A  =  E  in  the  first  integral  and  x  +  A  =  -E  in  the  second 
transforms  (A3. 2)  into 
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^NS 


(x  +  A)  :f[x  +  A  -  eV)-f  (x  +  A)|  dx 


[x  (x  +  2A)i 


',1/2 


+ 


(A3. 3) 


(x  +  Aj  I  f  (-X  -  A  -  eV)  -  f  (-  X  -  A)'  dx 
(x  +  2A)i’ 


\ll/2 


Now  using  the  general  relationship  valid  for  the  Fermi  function 


(A3.4) 


and  carrying  out  the  algebra  gives 


(x  +  A)  f(x  +  A- eV)-f  (x  +  A  +  eV)i  dx 
[x  (x  +  2A)]’^ 


(A3. 5) 


This  hasn't  quite  removed  the  singularity  but  allowing  the  substitution  of 
u2  =  X  will.  Rewriting  (A3. 5)  with  this  substitution  and  then  simply 
changing  the  dummy  variable  back  to  x  yields 


*NS 


2|x^  +  A|[f(x  +  A  -  eV)-f(x  +  A  +  eV)j  dx 


x^  +  2A 


il/2 


(A3. 6) 
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Finally,  this  integrand  is  in  a  convenient  form  with  which  to  complete  the 
numerical  integrations.  It  is  with  (A3. 6}  that  the  included  FORTRAN  program 
executed,  to  include  the  accompanying  derivatives. 

The  D01AHE  -  NAG  FORTRAN  library  routine  computes  a  definite  integral 


/, 


of  the  form  I  f(x)dx  .  The  method  uses  as  its  basis  a  family  of  interlacing 


high  precision  rules  as  outlined  in  reference  [41].  Two  advantages  of  this 
specific  subroutine  include; 

(a)  The  rate  of  convergence  is  monitored  and  used  to  make  a  decision  to 
abort  and  subdivide  before  the  full  sequence  has  been  applied. 

(b)  The  relative  accuracy  sought  in  each  subinterval  is  adjusted  in 
accordance  with  its  likely  contribution  to  the  total  integral. 

The  subroutine  parameters  and  call  are  as  follows: 


real  FUNCTION  D01AHE  (A,B,EPSR,NPTS,RELERR,F,NLIMIT,IFAIL) 
INTEGER  NPTS,  NLIMIT,  IFAIL 
REAL  A,B,EPSR,RELERR,F 
EXTERNAL  F 

where, 

A  —  lower  bound  of  integral 
B  --  upper  bound 

EPSR  "  specifies  relative  accuracy  required 
NPTS  -  counter  for  the  number  of  evaluations  used 
RELERR  contains  the  relative  error  achieved  (rough  guess) 

F  -  is  the  function  supplied  by  the  user 

NLIMIT  -  specifies  the  limit  of  the  number  of  function  evaluations: 
default  is  10,000 

IFAIL  --  is  a  flagged  fail  message 
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Program  FIT5 


c 

C  This  version  computes  I (V) ,  dliv,  d2iv,dliT  (correctly!)  all  with 
C  delta  a  function  of  temperature.  The  d2iT  calculations  have  been 
C  removed.  Tc  is  assigned  inside  the  program,  not  as  an  interactive 
C  input.  In  addition,  this  computes  d3idv3 . 

C 

Common/TJ/  delta,  T,  v,  beta,  SQtgap,  Tc 
Real  D2iVo(1000),  Volts(lOOO),  Predi(lOOO) 

Real  DliVo(lOOO),  DliTe(lOOO),  D3iVo(1000) 

Real  Cerror { 1000) ,  Verror(lOOO) 

Real  vstep,  vrangel,  vrange2,  EPSR,  RELERR,  high 
Real  DOIAHE 

Integer  IFAIL,  NLIMIT,  I 

External  current,  d2iv,  dliv,  dliT,  d3idv3 
C 
C 

c 

NLIMIT  =  0 
EPSR  =  1 .OE-6 
IFAIL  =  0 
C 

C2345678 

C 

Write{*,*)  'Input  T  (should  NOT  be  higher  than  Tc  =  6.79768) 
Read (*,10)  T 
10  Format (FIO .4) 

Write {*,*)  'Input  Vstep  size:  ' 

Read (*,10)  vstep 

Write(*,*)  'Beginning  Voltage(mV):  ' 

Read (*,10)  vrangel 

Write(*,*)  'Ending  voltage (mV) :  ' 

Read (*,10)  vrange2 
C 

C2345678 

C 

high  =  (vrange2  -  vrangel)  /  vstep 
beta  =  1 . /  (  . 086  *  T) 

Tc  =  6.79768 

SQtgap  =  SQRTd  -  (T/Tc)  ) 
delta  =  1.6  *  0.086  *  Tc  *  SQtgap 
C 
C 
C 

Open(unit  =  16,  file  =  ' curves . out ' ,  form  =  'formatted', 

$  status  =  'new') 

C  Open(unit  =  17,  file  =  ' errors . out ’ ,  form  =  'formatted', 

C  $  status  =  'new' ) 

C 

C 

c 

Write(*,20)  T,  delta,  vstep 
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c 

20 

$ 

30 

$ 

C 

C 

c 

c 


$ 


$ 


$ 

$ 

$ 

35 
C 
C 

C2345678 

Write(16,40)  (Voltsd),  Predid),  D2iVod), 

$  DliVod),  DliTed)  ,D3iVod)  ,  1  =  1,  high+1) 

40  Format(E13.7,  IX,  E13.7,  IX,  E13.7,  IX,  E13 . 7 , IX, E13 . 7 , 

$  lx,E13.7) 

C 

C 

C 

C  Writed7,50) 

C50  Format{ 'Volts (mV) ' ,1X, ' (C)error' ,1X, ' (V)error' ) 

C  Writed7,40)  (Volts  (I )  ,  Cerror  (I )  ,  Verrord), 

C  $  1=1,  high  +  1) 

C 

C 

Close  (unit  =  16,  status  =  'keep') 

C  Close  (unit  =  17,  status  =  'keep') 

C 

C 

End 

C 

C 

C2345678 

Real  Function  current (x) 

Common/TJ/  delta,  T,  v,  beta,  SQtgap.  Tc 
Real  x,gl,g2,sn 


Writed6,20)  T,  delta,  vstep 
Write (17, 20)  T,  delta,  vstep 

Format (' Temperature=  ' , FIO . 4/ , ' Delta=  ' , FIO .4/ , 'Volt  Step= 

F10.4//) 

Write (16,30) 

Format ( 'Volts (mV) ' , IX, 'Current' , IX, ' 2d-Deriv(V) • , iX, 

' Ist-Deriv(V) ' , IX, ' Ist-Deriv(T) ' , IX, ' 3rd-Deriv(V) ' ) 


DO  35  I  =  1,  high+1 
V  =  (1-1)  *  vstep  +  vrcingel 
Voltsd)  =  V 

Predid)  =  DOIAHE  (0 .  ,vrange2  +  3 ,  EPSR,  NPTS,  RELERR, 
cur rent, NLIMIT, IFAIL) 

Cerror (I)  =  RELERR 

D2iVod)  =  D01AHE(0  .  ,vrange2+3  .  ,  EPSR,NPTS,RELERR, 
d2iv,  NLIMIT, IFAIL) 

Verrord)  =  RELERR 

DliVod)  =  DO  lAHE  ( 0.  ,vrange2+3.  ,  EPSR,  NPTS,  RELERR, 
dliv,NLILMIT, IFAIL) 

DliTed)  =  DOIAHE  (  0.  ,vrange2+3.  ,  EPSR,  NPTS,  RELERR, 
dliT, NLIMIT, IFAIL) 

D3iVo(I)  =  D01AHE(0. ,vrange2+3. , EPSR, NPTS, RELERR, 
d3idv3, NLIMIT, IFAIL) 

CONTINUE 
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c 

c 

c 

c 


c 

c 

c 


gl  =  exp  (beta  *  ( (x**2)  +  delta  -  v)  ) 

g2  =  exp  (beta  *  ((x**2)  +  delta  +  v)  ) 

sn  =  ((x**2)  +  delta)  /  SQRT((x**2)  +  2*delta) 

current  =2*sn*  (  (1/(1+  gl)  )  -1/  (l+g2)  ) 

Return 

End 


Real  Function  d2iv(x) 

Commcn/TJ/  delta,  T,  v,  beta,  SQtgap,  Tc 

Real  V, delta, x,gl,g2.g3,g4,sn 

gl  =  exp  (beta  *  ( (x**2)  +  delta  -  v)  ) 

g2  =  exp  (beta  *  ((x**2)  +  delta  +  v)  ) 

sn  =  ((x**2)  +  delta)  /  SQRT((x**2)  +  2*delta) 

g3  =  gl  *  (gl-1)  /  (1  +  gl)**3 

g4  =  g2  *  (g2-l)  /  (1  +  g2)**3 

d2iv  =  2  *  sn  *  (beta**2)  *  (g3  -  g4) 

Return 

End 


Real  Function  dliv(x) 

Common /TJ/  delta,  T,  v,  beta,  SQtgap,  Tc 

Real  T,  V,  delta, x,gl,g2,sn 

gl  =  exp  (beta  *  ((x**2)  +  delta  -  v)  ) 

g2  =  exp  (beta  *  ( (x**2)  +  delta  +  v)  ) 

sn  =  ((x**2)  +  delta)  /  SQRT((x**2)  +  2*delta) 

dliv  =  2*sn*beta*  ( (gl/ (1+gl) **2)  +  g2/ (l+g2/ **2) 

Return 

End 


C 

C2345678 

C 

Real  Function  dliT(x) 

Common/TJ/  delta,  T,  v,  beta,  SQtgap,  Tc 

Real  T,  V,  delta, x,gl,g2,sn 

gl  =  exp  (beta  *  ( (x**2)  +  delta  -  v)  ) 

g2  =  exp  (beta  *  ( (x**2)  +  delta  +  v)  ) 

sn  =  ((x**2)  +  delta)  /  SQRT((x**2)  +  2*delta) 

C 

C  The  following  gobbledegoo)c  is  the  1st  derivative  of  current  w.r.t 
C  temperature . 

C 

C2345678 

dliT  =  -0.1376* (1/ (1  +  gl)  -  1/(1  +  g2))/ 

$  (SQtgap*Sqrt (0 . 2752*SQtgap*Tc  +  x**2))  + 

$  0.1376*(1/(1  +  gl)  -  1/(1  +  g2))*(x**2  +  delta)/ 

$  (SQtgap* (0.2752*SQtgap*Tc+  x**2) *• (3/2) )  + 

$  2*sn* (- (gl* (-0.8/ (T*SQtgap)  - 

$  11.6279*(-v  +  x**2  +  delta) /T**2) / (1  +  gl)**2) 

$  +  g2* (-0.8/ (T*SQtgap)  - 
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$  11.6279*(v  +  x**2  +  delta) /T**2) / (1  +  g2)**2) 

Return 

End 

c 

C2345678 

C 

Real  Function  d3idv3 (x) 

Common/TJ/  delta,  T,  v,  beta,  SQtgap,  Tc 
Real  T,  V,  delta,x,subl,sub2,TK3,sn 

subl  =  exp ( (1 . 6*0 . 086*SQtgap*Tc  +  v  +  x**2)  /  (.086*T)) 
sub2  =  exp ( (1.6*0. 086*SQtgap*Tc  -  v  +  x*2)  /  (.086*T)) 
sn  =  ((x**2)  +  delta)  /  SQRT((x**2)  +  2*delta) 

TK3  =  (0.086*T)**3 
C 

C  The  following  gobbledegoo)c  is  the  3rd  derivative  of  current  w.r.t 
C  voltage. 

C 

C23456789 

d3idv3  =  2*sn* ( (6*subl**3) / ( (1  +  subl) **4*TK3)  - 
$  6*subl**2/ ( (1  +  subl) **3*TK3)  + 

$  subl/{(l  +  subl) **2*TK3)  + 

$  (6*sub2**3) /  I  (1  +  sub2) **4*TK3)  - 

$  6* (sub2**2) / ( (1  +  sub2) **3*TK3)  + 

$  sub2/({l  +  sub2) **2*TK3) ) 

Return 

End 


86 


References 


1  G.M.  Elchinger,  "Properties  of  Metal-Barrier-Metal  Tunneling 
Junctions  Responding  to  Infrared  and  Visible  Radiations",  (PhD. 
thesis,  Mass.  Inst,  of  Tech.,  1978) 

2  K.C.  Liu,  "Characteristics  of  Point  Contact  Submicron  Junctions  and 
Their  Applications",  (PhD.  thesis.  Mass.  Inst,  of  Tech.,  1979) 

3  L.  Solymar,  Superconductive  Tunnelling  and  Applications,  (Wiley- 
Interscience,  New  York,  1972) 

4  L.O.  Mocker,  and  A.  Javan,  "Laser  Harmonic  Frequency  Mixing  of  Two 
Different  Far  Infrared  Laser  Lines  up  to  118pm",  Phys.  Lett.  26 A, 
255-6  (1968) 

5  L.O.  Mocker,  D.R.  Sokoloff,  V.  Daneu,  A.  Szoke,  and  A.  Javan, 

"Frequency  Mixing  in  the  Infrared  and  Far-Infrared  Using  a  Metal-to- 
Metal  Point  Contact  Diode",  App.  Phys.  Lett.  12,  401-2  (1968) 

6  M.  Heiblum,  S.Y.  Wang,  J.R.  Whinnery,  and  T.K.  Gustafson, 
"Characteristics  of  Integrated  MOM  Junctions  From  DC  to  Optical 
Frequencies",  J.  Quantum  Elect.  13,  48  (1977) 

7  J.G.  Small,  G.M.  Elchinger,  A.  Javan,  A.  Sanchez,  F.J.  Bachner,  and  D.L. 
Smythe,  "AC  Electron  Tunneling  at  Infrared  Frequencies:  Thin  Film  M- 
0-M  diode  Structures  With  Broad-Band  Characteristics",  App.  Phys. 
Lett.  24,  275  (1974) 

8  T.K.  Gustafson,  R.V.  Schmidt,  and  J.R.  Perucca,  "Optical  Detection  in 
Thin  Film  Metal-Oxide-Metal  Diodes",  App.  Phys.  Lett.  24,  620  (1974) 

9  S.Y.  Wang,  T.  Izawa,  and  T.K.  Gustafson,  "Coupling  Characteristics  of 
Thin-Film  Metal-Oxide-Metal  Diodes  at  10.6  pm",  App.  Phys.  Lett. 
27,481  (1975) 

10  B.  Fan,  S.M.  Farris,  T.K.  Gustafson,  and  T.J.  Bridges,  "Nonlinear  Optical 
Response  of  Metal-Barrier-Metal  Junctions",  App.  Phys.  Lett.  30,  177 
(1977) 


87 


1 1  Richard  L.  Liboff,  Introductory  Quantum  Mechanics,  (Addison-Wesley 
Publishing  Company,  New  York, 1989) 

12  Claude  Cohen-Tanoudji,  Bernard  Diu,  and  Franck  Laloe,  Quantum 
Mechanics,  (John  Wiley  and  sons.  New  York,  1977) 

13  J.G.  Simmons,  "Generalized  Formula  for  the  Electric  Tunnel  Effect 
Between  Similar  Electrodes  Separated  by  a  Thin  Insulating  Film",  J. 
App.  Phys.  34,  1793-1803  (1963) 

14  R.  Stratton,  "Volt-Current  Characteristics  for  Tunneling  Through 
Insulating  Films",  J.  Phys.  Chem.  Solids  23,  1177-1190,  (1962) 

1 5  E.O.  Kane,  "  Basic  Concepts  of  Tunneling",  Tunneling  Phenomena  in 
Solids,  E.  Burstein  and  S.  Lundquist,  eds.,  (Plenum  Press,  New  York, 
1969) 

16  J.  Bardeen,  Phys.  Rev.  Letts.  6,  57-9  (1961) 

17  I.  Giaever,  "Energy  Gap  in  Superconductors  Measured  by  Electron 
Tunneling",  Phys.  Rev.  Lett.  5,  147-8  (1960) 

18  J.C.  Fischer  and  I.  Giaever,  "Tunneling  Through  Thin  Insulating 
Layers",  J.  App.  Phys.  32,  172-177  (1961) 

19  I.  Giaever  and  K,  Megerle,  "Study  of  Superconductors  by  Electron 
Tunneling",  Phys.  Rev.  Letts.  14,  108  (1965) 

20  A.  Dayem  and  R.J.  Martin,  "Quantum  Interaction  of  Microwave 
Radiation  with  Tunneling  Between  Superconductors",  Phys.  Rev. 
Letts.  8,  246  (1962) 

21  P.K.  Tien  and  J.  Gordon,  "Multiphoton  Process  Observed  in  the 
Interaction  of  Microwave  Fields  with  the  Tunneling  Between 
Superconductor  Films",  Phys.  Rev.  129,  647  (1963) 

22  Y.  Goldstein  and  B.  Abeles,  "Superconducting  Tunneling  Induced  by 
Gigacycle  Sound  Waves",  Phys.  Letts.  14,  78  (1965) 

23  E.  Lax  and  F.L.  Vernon,  Jr.,  "Microwave  Phono-Assisted  Tunneling  in 
Superconductors",  Phys.  Rev.  Letts.  14,  256  (1965) 

24  B.N.  Taylor  and  E.  Burstein,  "Excess  Currents  in  Electron  Tunneling 
Between  Superconductors",  Phys.  Rev.  Letts.  10,  14  (1963) 


88 


25  C.J.  Adkins,  "Two  particle  Tunneling  Between  Superconductors",  Phil. 
Mag.  8,  1051  (1963) 

26  J.M.  Rowell,  "Tunneling  Between  Superconductors",  Rev.  Mod.  Phys. 

36,  215  (1964) 

27  I.K.  Yanson,  V.M.  Svistunov  and  I.M.  Dmitrenko,  "Experimental 
Observations  of  Tunneling  of  Cooper  Pairs  Between  Thin  Layers  of 
Superconducting  Tin",  JETP  20,  1404  (1965) 

28  W.J.  Tomasch,  "Geometrical  Resonance  in  the  Tunneling 
Characteristics  of  Superconducting  Pb",  Phys.  Rev.  Letts.  15,  672 
(1965) 

29  Z.  Burshtein  and  J.  Levinson,  "Photo-induced  Tunnel  Currents  in  Al- 
AI2Q3-AU  Structures",  Phys.  Rev.  B  12,  3453-7  (1975) 

30  R.C.  Jaklevic  and  J.  Lambe,  "Molecular  Vibration  Spectra  By  Electron 
Tunneling",  Phys.  Rev.  Lett.  17,  1139-40  (1966) 

31  A.  Braunstein,  M.  Braunstein,  and  G.S.  Picus,  "  Photoemissive 
Determination  of  Barrier  Shape  in  Tunnel  Junctions",  Phys.  Rev.  Lett. 
14,  219-21  (1965) 

32  G.M.  Elchinger,  A.  Sanchez,  C.F.  Davis,  and  A.  Javan,  "Mechanism  of 
Detection  of  Radiation  in  a  High-Speed  Metal-Metal  Oxide-Junction 
in  the  Visible  Region  and  at  Longer  Wavelengths",  J.  App.  Phys.  47, 
591-4  (1976) 

33  A.  Sanchez,  C.F.  Davis,  G.  Lui,  and  A.  Javan,  "The  MOM  Tunneling  Diode; 
Theoretical  Estimate  of  its  Performance  at  Microwave  and  Infrared 
Frequencies,"  J.  AppI  Phys.  49,  10  (1978) 

34  J.  Vrba,  and  S.B.  Woods,  "Glow-Discharge  of  Freshly  Evaporated 
Aluminum  Films",  Can.  J.  Phys.  50,  548-52  (1972) 

35  J.L.  Miles,  and  P.H.  Smith,  "The  Formation  of  Metal  Oxide  Films  Using 
Gaseous  and  Solid  Electrolytes",  J.  Electrochem.  Soc.  110,  1240-5 
(1963) 

36  D.E.  Thomas  and  J.M.  Rowell,  "Low  Level  Second  Harmonic  Detection 
System,"  Rev.  Sci.  Instr.  36,  9  (1965) 

37  M.H.  Cohen,  L.M.  Falicov  and  J.C.  Phillips,  "Superconductive 
Tunneling",  Phys.  Rev.  Letts.  8,  316  (1962) 


89 


J.R.  Schrieffer,  "  Single  Particle  Tunneling  in  Superconductors". 
Tunneling  Phenomena  in  Solids,  E.  Burstein  and  S.  Lundquist,  eds., 
(Plenum  Press,  New  York,  1969) 

J.R.  Schrieffer,  "Theory  of  Electron  TunnelingRev.  Mod.  Phys.  36,  200 
(1964) 

J.  Bardeen,  L.N.  Cooper,  and  J.R.  Schrieffer,  "Theory  of 
Superconductivity",  Phys.  Rev.  108,  5,  1175-1204  (1957) 

T.N.L.  Patterson,  "The  Optimum  Addition  of  Points  to  Quadratrue 
Formulae",  Math.  Comp.  22,  847-856  (1968) 

K. W.  Sheppard,  "Photocurrents  Through  Thin  Films  of  AI2O5'"  J.  App. 
Phys.  36,  796  (1965) 

J.D.  Jackson,  Classical  Electrodynamics,  (John  Wiley  and  Sons,  New 
York,  1975) 

V.E.  Griineisen,  "Die  Abhangigkeit  des  elektrischen  Widerstandes 
reiner  Metalle  von  der  Temperatur",  Ann.  d.  Phys.  5,  16,  530-540 
(1933) 

H.A,  Boorse,  and  H.  Niewodniczanski,  "The  Electrical  Resistance  of 
Aluminum  at  Low  Temperatures",  Proc.  Roy.  Soc.  A  153,  463  (1936) 

J. G.  Thomas  and  E.  Mendoza,  "LXXXVI.  The  Electrical  Resisance  of 
Magnesium,  Aluminum,  Molybdenum,  Cobalt,  and  Tungsten  at  Low 
Temperatures",  Phil.  Mag.,  900-910  (1951) 

D.  Chaverri,  A.  Saenz,  and  V.  Castano,  "Grain  Size  and  Electrical 
Resistivity  Measurements  on  Aluminum  Polycrystalline  Thin  Films", 
Materials  Letters  12,  344-348  (1991) 

K. M.  Evenson,  D.A.  Jennings,  F.R.  Peterson,  and  J.S.  Wells,  Laser 
Spectroscopy  III,  J.L.  Hall  and  L.L.  Carlsten,  eds.,  (Springer-Verlag, 
Berlin,  1977).  pp  56-68. 

L. O.  Mocker,  J.G.  Small,  and  A.Javan,  "Extension  of  Absolute  Frequency 
Measurements  to  the  84p  Range",  Phys.  Letts.  29 A,  6,  321-2  (1969) 

E. N  Economou,  "Surface  Plasmons  in  Thin  Films",  Phys.  Rev.  182,  2, 
539  (1969) 


